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ABSTRACT 


This  research  considers  the  devising  of  climatic 
presentations  for  usage  in  preparing  2-  and  4-hour 
forecasts  of  low  ceilings  and  restricted  visibilities 
and  the  providing  of  integrated-wind  information  within 
the  drop-zone  layer  in  support  of  paradrop  operations. 

A  hierarchy  of  empirical  relationships  for  esti¬ 
mating  the  recurrence  probabilities  under  a  number  of 
contingencies  are  made  possible  through  the  use  of 
statistical  modelling  techniques  which  act  to  conserve 
the  data  base.  These  are:  1)  a  model  for  generating 
recurrence  probabilities  which  considers  only  the  latest 
observed  ceiling  category,  2)  one  which  incorporates 
surface  wind  information,  3)  another  which  includes 
moisture  considerations  and  4)  one  which  introduces  ob¬ 
served  moisture  and  wind  information  diagnostically  with¬ 
out  having  to  assume  the  persistency  of  these  variables, 
or  otherwise  forecast  their  behavior,  throughout  the 
forecast  interval. 

Climatic  presentations  for  estimating  integrated 
winds  within  the  drop-zone  layer,  from  a  knowledge  of  the 
winds  at  flight  level,  are  discussed  apart  from  those 
which  pertain  when  surface  winds  are  similarly  available. 
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I.  INTRODUCTION 


Early  in  the  program  the  AFCRL-AWS-SLU  team  com¬ 
pleted  the  first  task  of  the  contract—  to  investigate 
the  critical  problem  areas  to  be  researched  and  estab¬ 
lish  criteria  and  priorities  for  attacking  them.  Two 
specific  areas  were  outlined:  1)  Devising  climatic 
presentations  for  operational  Air  Weather  Service  (AWS) 
usage  in  preparing  2-  and  4-  hour  forecasts  of  low 
ceilings  and  restricted  visibilities,  and  2)  designing 
climatic  information  to  support  paradrop  operations 
under  conditions  when  only  the  winds  at  flight  and 
ground  levels  are  known. 

The  ground  rules  for  research  developments  were 
also  established  at  that  time.  First  priority  was  to 
be  given  to  the  undertaking  of  feasibility  studies  to 
establish  the  relative  merit  of  each  proposed  attack. 

Once  feasibility  was  established,  solutions  were  to  be 
sought  which  offer  possibilities  of  generalized  appli¬ 
cability  throughout  the  diverse  areas  of  AWS  interests. 
The  final  products  were  to  be  presented  in  formats 
acceptable  to  the  operational  AWS  forecaster.  All  the 
problems  were  to  be  considered  from  a  climatological 
point  of  view,  i.e.,  the  contract  calls  for  providing 
climatic  information  to  assist  the  forecaster — not 
designing  it  to  eliminate  him.  The  research  efforts  for 
the  two  aforementioned  areas  will  be  presented  sepa¬ 
rately  in  this  report. 

The  problem  of  forecasting  low  ceilings  is  not  new 
to  meteorology.  In  fact,  the  literature  abounds  with 
reports  of  partial  successes  and  failures  by  others  who 
have  attempted  to  solve  the  problem.  Hence,  our  first 
task  was  to  review  these  studies  in  order  to  circumvent 
as  many  "false  starts"  as  possible  during  the  course  of 
our  research.  It  is  perhaps  trivial  to  state  that  the 
problem  would  have  been  solved  long  before  if  ready 
solutions  existed.  However,  this  stark  realization  is 
important  from  the  outset  in  order  to  avoid  discourage¬ 
ments  later.  Consequently,  the  research  was  oriented 
toward  the  devising  of  climatic  presentations  to  be 
used  as  forecast  aids  rather  than  the  seeking  of  fore¬ 
cast  entities  which,  in  all  probability,  are  beyond  our 
abilities  to  produce  for  a  number  of  reasons.  Foremost 
among  these  is  simply  the  fact  that  the  data  records  are 
not  sufficiently  refined  to  contain  many  of  the  short 
term  meteorological  fluctuations  which  affect  the  ceiling 
conditions  over  this  time  period.  The  data  problem  is 
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particularly  significant  for  low  ceiling  heights  due  to 
the  added  complications  imposed  by  the  local  geography 
and  topography  of  the  earth's  surface.  For  example,  on 
occasions  low  ceilings  are  produced  by  local  effects  and 
on  the  others  by  quasi- stationary  weather  systems  which 
cover  wide  areas  of  the  map.  Low  ceilings  may  also  occur 
in  conjunction  with  well-defined  migratory  features. 

Thus,  at  times  the  best  procedures  for  forecasting  ceiling 
changes  may  be  extrapolative  or  advective  in  character 
while  on  others  the  best  tools  to  use  may  be  simple  per¬ 
sistency  or  its  climatology. 

AWS  has  devised  persistence  probability  tables  (PPT) 
which  are  particularly  applicable  under  stagnant  weather 
regimes.  We  chose  to  build  upon  this  expertise  and  use 
their  processed  data  bases  for  devising  and  testing  im¬ 
proved  applications  of  this  basic  concept. 

Since  climatology  is  but  one  approach  to  the  problem 
of  forecasting  ceilings,  the  forecaster  is  encouraged  to 
use  its  guidance  along  with  that  provided  by  a  synthesis 
of  the  other  materials  at  his  disposal.  The  relative 
weighing  of  these  forecast  informations  is  a  judgment 
factor  which  he  is  best  suited  to  make  since  it  depends 
upon  the  meteorological  conditions  under  which  a  particular 
ceiling  is  occurring  and  their  past  histories  of  movement 
and  association.  Our  research  aim  is  to  incorporate  ob¬ 
jectively  the  more  common  parameters  which  affect  ceiling 
persistency  so  that  the  problem  becomes  primarily  one  of 
deciding  when  to  use  climatology  rather  than  how  to  use  it. 

Most  researchers  who  have  attacked  this  problem  have 
relied  heavily  on  statistical  techniques  to  uncover  fore¬ 
cast  relationships  of  note.  We  chose  to  do  likewise  and 
solicit  the  aid  of  the  computer  to  diagnose  the  physical 
significance  of  these  relationships.  In  particular,  we 
hope  to  incorporate  objective  considerations  based  upon 
time  of  day,  season  of  the  year,  wind  speed  and  direction, 
and  humidity  conditions  into  the  climatic  recurrence 
statistics . 

A  most  important  consideration  which  must  be  made  at 
the  outset  is  how  best  to  display  these  statistical  re¬ 
lationships.  If  the  approach  is  to  involve  a  mix  between 
the  data  processing  attributes  of  the  computer  and  the 
interpretative  abilities  of  man,  these  displays  must  nec¬ 
essarily  be  kept  sufficiently  simple  to  permit  an  effective 
interplay.  Many  of  the  failures  in  the  past  have  resulted 
from  inabilities  to  recognize  that  man  has  this  role  to 
play.  Either  by  direct  design,  or  because  the  statistical 
methods  being  employed  rapidly  became  too  involved  for  man 
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to  mentally  interpret,  the  forecaster  has  been  handed 
categorical  forecasts  produced  by  a  combination  of  re¬ 
gression  terms  which  he  is  asked  to  either  accept  or 
reject,  in  total,  with  little  guidance  on  when  to  do 
either.  To  avoid  this  pitfall,  we  chose  to  use  simple 
diagrammatic  techniques  which  pictorially  portray 
various  non-linear  combinations  of  meteorological  param¬ 
eters  under  synoptic  conditions  which  are  most  commonly 
found  in  the  natural  atmosphere. 

Having  decided  on  the  problem  to  attack  and  the 
statistical  approaches  to  use,  we  turned  our  attention 
to  the  most  formidable  problem  yet  to  be  resolved.  How 
can  the  problem  of  ceiling  persistency  best  be  resolved 
in  view  of  the  restrictions  that  a  limited  data  base 
imposes  upon  it?  In  attempting  to  answer  this  question, 
we  have  systematically  progressed  through  four  complexi¬ 
ties  of  statistical  modelling.  These  are:  1)  a  model 
for  generating  recurrence  probabilities  which  considers 
no  meteorological  input  other  than  a  knowledge  of  the  ob¬ 
served  ceiling  category  at  the  time  of  forecast,  2)  a 
model  which  incorporates  the  latest  observed  surface 
wind  information,  3)  a  version  which  includes  moisture  con¬ 
siderations  (the  latest  observed  temperature  dew-point 
spread),  and  4)  a  revised  concept  which  incorporates 
moisture  and  wind  information  diagnostically  without  re¬ 
quiring  that  any  prognostic  information  whatsoever  be 
supplied . 

II.  GENERATING  RECURRENCE  PROBABILITIES  BY  A  MODEL  WHICH 
CONSIDERS  NO  METEOROLOGICAL  INPUT  AT  THE  TIME  OF 
FORECAST 

Probabilities  that  given  ceiling  events  will  recur 
at  later  times  can  be  obtained  by  various  means.  Those 
used  as  a  standard  for  comparisons  and  for  model  de¬ 
velopments  in  this  research  were  provided  by  AWS.  They 
were  computer-generated  by  that  establishment  in  its 
Asheville,  North  Carolina  facility  by  simply  tracing  each 
event  throughout  the  hourly  records.  These  were  made 
available  to  us  for  a  number  of  ceiling  categories  and 
stations  located  throughout  the  world.  They  provide  a 
direct  evaluation  of  persistency  without  resorting  to  any 
assumptions  whatsoever.  They  are  as  reliable  as  the  data 
records  will  provide  by  "brute -force"  computer  processings. 
Therefore,  if  sufficient  data  were  available  to  insure 
statistical  significance  to  the  results,  and,  if  computer 
processing  time  did  not  constitute  a  restraining  influence, 
there  would  be  no  need  to  resort  to  the  modelling  pro¬ 
cedures  as  we  are  about  to  outline.  A  fundamental  problem 
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arises.  Sufficient  data  are  not  available  anywhere,  for 
any  station,  to  incorporate  objectively  the  many  meteoro¬ 
logical  considerations  which  the  forecaster  is  presently 
attempting  to  do  subjectively  in  his  assessments  of  the 
persistency  qualities  of  a  given  observed  event. 

The  model  offers  a  means  of  circumventing  some  of 
those  processing  steps  which  rapidly  consume  the  data 
base  by  relying  on  certain  other  well-founded  statistical 
relationships.  McCabe  (1968)  estimates  that  it  requires 
ten  times  as  much  data  to  produce  reliable  conditional 
probabilities  that  an  event  will  reoccur  as  it  does  to 
determine  their  unconditional  frequency  of  occurring  in 
the  first  place.  Hence,  the  specific  step  which  we  hope 
to  circumvent  involves  the  necessity  of  computer  tracing 
the  subsequent  behavior  of  each  given  observed  event. 

Only  by  conserving  data  in  this  way  are  we  able  to  have 
enough  left  to  incorporate  consideration  of  the  embedding 
synoptic  environment  in  which  a  given  ceiling  is  being 
observed  into  our  recurrence  probability  statistics. 

Three  specific  models  were  found  in  the  literature  which 
had  a  direct  bearing  on  our  problem.  They  were  provided 
by  McCabe  (1968),  McAllister  (1969),  and  Gringorten  (1971). 
We  chose  to  evaluate  the  latter  two  prior  to  embarking  on 
the  developing  of  models  of  our  own.  These  two  respective 
evaluations  appear  to  be  a  relatively  simple  matter.  Each 
uses  a  pair  of  unconditional  frequency  of  ceiling  occur¬ 
rence  climatologies  (one  for  to  and  the  other  for  tQ  +  2 
in  the  case  of  a  2-hour  forecast)  as  input  parameters  to 
generate  conditional  climatologies  of  recurrence  proba¬ 
bilities.  Thus,  both  the  necessary  input  and  output 
values  for  these  evaluations  are  contained  in  the  AWS  com¬ 
puter-generated  data.  However,  there  is  one  big  problem. 
The  raw  statistics  provided  by  AWS  are  not  consistently 
reliable  for  the  simple  reason  that  during  certain  times 
of  the  day  and  months  of  the  year  an  insufficient  number 
of  low-ceiling  occurrences  were  observed  to  provide  a 
reliable  set  of  recurrence  statistics  by  the  direct 
computer-tracing  approach.  Also,  ceilings  are  often 
rendered  fictitiously  discontinuous  (by  the  very  defini¬ 
tion  of  a  ceiling)  since  wild  fluctuations  in  ceiling 
heights  are  reported  whenever  the  accumulated  cloud  cover- 
ate  oscillates  between  broken  and  scattered.  Fortunately 
many  of  these  effects  can  be  minimized  by  subjecting  the 
climatic  representations,  which  are  derived  from  these 
individual  data  reports ,  to  hour-to-hour  and  month-to- 
month  time  continuity  constraints.  The  need  for  some  such 
procedure  for  "cleaning-up"  the  data,  without  destroying 
its  information  content,  becomes  increasingly  apparent 
whenever  one  shrinks  the  original  data  base  by  stratifying 
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it  into  "homogeneous"  subsets  of  meteorological  occur¬ 
rences.  The  validity  of  using  some  such  continuity 
constraints  was  confirmed  by  the  data  themselves  when¬ 
ever  the  probability  statistics  were  derived  from  suf¬ 
ficiently  large  samples  to  dwarf  the  errors  arising 
from  the  aforementioned  more  random  considerations.  In¬ 
variably  we  found  that  the  more  frequent  the  occurrence 
the  better  the  data  fit  the  climatic  time-continuity 
constraints.  The  format  of  Figure  1  was  designed  to  en¬ 
hance  and  maintain  the  integrity  of  the  data  as  a  repre¬ 
sentative  climatology  by  these  considerations. 

hour  n  ?  *  *  *  in  12  1*  16  1«  20  ?2 


Fig.  1.  Observed  unconditional  frequency  of  ceiling 

occurrences  equal  to  or  less  than  200  feet  for 
Randolph  AFB,  Texas. 

Values  aligned  in  the  format  of  Figure  1  which  do 
not  conform  to  a  smoothed  isoline  of  frequency  of  occur¬ 
rence  probabilities  are  immediately  isolated  for  man's 
interpretation  by  providing  him  a  means  of  viewing  the 
misfits  against  the  backdrop  of  a  "continuous"  field  of 
meteorological  information.  The  field  analysis  of  a 
station's  climatic  records  for  a  given  set  of  occurrence 
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frequencies  can  be  best  performed  by  systematic  pro¬ 
gressions  from  the  more  reliable  (highest  frequency  of 
occurrence)  statistics  toward  those  areas  occupied  by 
the  more  suspect  ones.  Meteorologists  are  particularly 
adept  in  employing  such  procedures  since  most  of  the 
analyses  they  have  been  trained  to  make  involve  mixes  of 
reliable  and  unreliable  data  over  dense  and  sparse 
coverages.  The  isoline  of  probability  analyses  can  either 
be  made  subjectively  or  objectively.  After  considerable 
experimentations  an  objective  3x3  smoothing  matrix 
(see  Figure  2)  was  found  to  give  values  at  the  respective 
grid  points  which  fit  those  produced  by  our  best  subjec¬ 
tive  analysts  within  accepted  tolerances.  The  data  in 
Figure  1  has  already  been  smoothed  in  the  above  manner. 

Fig.  2.  The  objective  2-d imeniional 
smoothing  operator 


Restriction:  The  statistics  for  any  event  must  change  smoothly 


and  continuously  across  the  2 
the  data. 


EXAMPLE: 


Actual  Data  June  12  LST 
699  Total  Observations 
17  Ceilings  <200  Feet  Month 
Probability  of 
Occurrence  =  .024 


dimensional  field  used  to  present 


Probability 
o F 

Ceiling  *200  ft. 


Time  of  Day 

To  smooth  data  for  any  month  and  time  (eg.  June  12  LST)  the 
0  surrounding  data  sets  are 


11  LST _ 12  LST _ 13  LST 


ttl - 

8T2 

— rr? - 

(20) 

(13) 

(14) 

.023 

.015 

.013 

1 

A 

1 

699 

699 

700 

(17) 

(17) 

(11) 

.  024 

.024 

.016 

B 

A.  B 

B 

860 

858 

859 

(24) 

(16) 

(14) 

.028 

.019 

.016 

1 

A 

1 

A  and  B  are  weighting  factors 

A  s  2  for  data  each  hour 
B  *  3  by  month  or  dew-point  spread 

B  =  3  if  dew-point  spread  replaces  month  on  vertical  axis 
SMOOTHED  UNCONDITIONAL  PROBABILITY  = 

(1)  (062)(  .023)  ♦  (2X862)  (.  015)  ♦  (  1)  (  .  862  )  ( .  013  >“J 
(3X699)  (  .024)  ♦  (  6  )  ( 6  99  )  (  .  024  )  ♦  (  3  )  (  .  7  00  )  (  .  01 6  ) 

( 1)  (  8  60 )  (  .  0  28  )  ♦  (2X858X  .019)  ♦  (1)(  .  8S0)(  .  016)  J 
( 1)  ( 8  6  2  )  ♦  2(  862  )  ♦  (1)  862)  ♦  (3)(699)  ♦  (6X699)  ♦  (3)(?',;7| 

(1X  860)  ♦  (  2X8  58  )  ♦  (1X  859)  J 

*  .0204 

Note  in  Figure  1  that  the  frequencies  of  low  ceiling  oc¬ 
currences  are  highly  dependent  on  the  local  hour-of-the 
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day  and  the  season-of-the  year  and  that  the  time  of  maximum 
and  minimum  occurrences  are  readily  perceived .  Particu¬ 
larly  noteworthy  is  the  fact  that  the  frequency  of 
ceiling  occurrences  is  not  predicated  on  any  forecast 
input  or  observed  knowledge.  Hence,  these  statistics  are 
presented  in  a  compact  format  that  should  be  particularly 
helpful  to  any  forecaster  who  assumes  duties  in  an  area 
for  the  first  time.  He  is  immediately  and  straight¬ 
forwardly  introduced  to  the  times  and  conditions  when  he 
is  most  likely  to  encounter  a  specific  forecast  problem 
at  that  location. 

Just  as  the  "observed"  unconditional  frequencies  of 
occurrences  were  found  to  change  smoothly  and  continu¬ 
ously  across  the  2-dimensional  display  of  Figure  1  so  did 
the  respective  "observed"  conditional  probabilities  which 
we  extracted  from  the  persistence  probability  tables  of 
AWS  (see  Figure  3).  This  display  discloses  at  a  glance 
the  hours  and  months  when  persistency  does  and  does  not 
offer  the  best  statistically  probable  forecast. 


Fig.  3.  Observed  2-hour  recurrence  probabilities  for 
.  ceilings  equal  to  or  less  than  200  feet  for 

Randolph  AFB,  Texas. 
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Formats ,  such  as  those  shown  in  Figures  1  and  3,  repre¬ 
sent  an  important  spin-off  of  this  research.  We  recom¬ 
mend  them  for  AWS  consideration  as  a  better  means  of 
presenting  reliable  climatic  information  than  the  ones 
now  in  being. 

The  information  contained  in  Figures  1  and  3 
readily  provides  a  sufficient  set  of  reliable  statis¬ 
tics  for  verifying  the  techniques  of  McAllister  and 
Gringorten.  McAllister  (1969)  investigated  the  recur¬ 
rence  probabilities  of  cirrus  clouds  and  found  that  "A 
sensitivity  analysis  shows  results  that  are  essentially 
independent  of  particular  data  source  with  regard  to 
geographic  location,  climatic  regime,  or  season  of  the 
year  insofar  as  parameter  determination  and  concomitant 
accuracy  of  estimation  are  concerned.  Extensions  to 
broader  classes  of  meteorological  phenomena  are  also  in¬ 
dicated."  His  findings  are  particularly  intriguing  since 
they  meet  one  of  our  aforementioned  objectives — to  orient 
the  research  towards  those  solutions  which  offer  the 
possibility  of  generalized  application  throughout  the 
areas  of  Air  Weather  Service's  interests.  An  evaluation 
of  McAllister's  model  was  undertaken  by  Captain 
Richard  Nieman,  (1971),  a  graduate  student  at  Saint  Louis 
University.  A  brief  discussion  of  his  research  follows: 

"The  empirical  recurrence  formula  proposed  by 
McAllister  was  used  to  estimate  unconditional 
probabilities  of  various  ceiling  categories 
for  three  stations  located  in  Texas.  These 
stations  lie  in  an  area  far  removed  and  climati¬ 
cally  different  from  the  area  that  McAllister 
used  to  derive  the  formula.  The  parameter 
values  given  by  McAllister  and  those  derived 
from  the  Texas  data  were  tested,  in  turn,  to 
determine  the  reliability  of  the  estimated 
conditional  probabilities  of  recurrence  under 
various  time,  area  and  ceiling  groupings.  The 
results  indicate  reliability  for  the  formula's 
estimates  under  conditions  of  large  unconditional 
probabilities.  The  method  breaks  down  when  the 
unconditional  probabilities  are  small .  Since  an 
exponentially  decreasing  function  is  involved, 
the  greatest  errors  in  the  estimates  occur  in  the 
shorter  forecast  periods.  Hence,  the  usefulness 
of  the  formula  is  dependent  on  the  frequency  of 
the  weather  phenomena  under  consideration  and  the 
time  period  of  the  forecast. 

Since  the  frequency  and  persistency  of  weather 
phenomena  are  dependent  on  factors  such  as  season, 
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time  of  day,  and  topography  in  combination 
with  dynamic  and  thermodynamic  processes, 
some  type  of  data  stratification  would  appear 
necessary  to  obtain  'unconditional  climatolo¬ 
gies'  which  are  less  biased  toward  those 
categories  which  are  highly  recurrent." 

Nieman's  conclusions  that  the  accuracy  of  McAllister's 
model  is  dependent  on  the  magnitude  of  the  unconditional 
probabilities  were  also  verified  by  Gringorten  (1971). 

These  investigations  show  that  McAllister's  model  is  not 
directly  applicable  to  our  problem  since  low  ceilings 
and  restricted  visibilities  occur  rather  infrequently  at 
most  stations  during  certain  hours  of  the  days  and  months 
of  the  year.  Nieman's  comments  on  stratification  were 
further  amplified  by  Major  Jerome  Wacker  (1971) ,  a  con¬ 
temporary  graduate  student  conducting  research  on  a 
allied  aspect  of  the  problem.  Wacker  showed  that  the  oc¬ 
currence  statistics  for  low  ceilings  at  Randolph  AFB, 

Texas  were  highly  dependent  on  the  wind  direction.  He 
concluded  that  data  stratification  effectively  "provides 
information  as  to  which  parameters  need  to  be  forecast, 
and  in  what  detail,  for  any  particular  station  or  problem 
area  where  a  solution  exists  within  the  available  cli¬ 
matic  records." 

A  model  set  forth  by  Gringorten  (1971)  was  also 
evaluated  using  the  aforementioned  computer-traced  statis¬ 
tics  of  AWS.  His  model  used  the  same  forms  of  input  data 
to  generate  conditional  probabilities  as  those  employed 
by  McAllister.  Gringorten  takes  advantage  of  the  Markov 
chain  assumption  and  the  U.  S.  Bureau  of  Standard's  bi¬ 
variate-normal  tables  to  produce  the  desired  result.  The 
unknown  variable  in  this  approach  is  the  hour-to-hour 
correlation  of  the  normalized  frequency  of  ceiling  occur¬ 
rences.  Its  value  is  assumed  in  his  model  to  be  0.95. 
Although  he  found  the  hour-to-hour  correlation  to  be  rela¬ 
tively  invariant  for  the  elements  which  he  tested, 
Gringorten  recognized  the  need  for  it  "to  be  surveyed  and 
mapped  as  a  function  of  element,  geography  and  possibility 
time  of  year" 

Data  of  the  type  presented  in  Figures  1  and  3  are 
sufficient  to  calculate  the  hour-to-hour  correlation 
value  explicitly  and  evaluate  its  constancy.  The  process 
merely  involved  the  input  of  the  known  2-hour  conditional 
probability  values  for  ceilings  less  than  200  feet  provided 
by  AWS  into  Gringorten' s  formula.  Such  calculations  based 
on  the  data  of  Randolph  AFB,  Texas  are  shown  in  Figure  4. 
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.959 
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.927 
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Fig.  4.  Hour  to  hour  correlation  values  ( f> )  for  ceilings 
less  than  200  ft.  at  Randolph  AFB,  Texas.  Calcu¬ 
late  from  Figures  1  and  3  and  bivariate  normal 
statistics . 


Here  the  values  are  found  to  be  hourly  and  seasonally 
variable.  The  question  arises,  Since  McCabe  (1968), 
McAllister  (1969)  and  Gringorten  (1971)  found  certain 
relationships  between  unconditional  and  conditional 
probabilities  which  were  relatively  independent  of  geo¬ 
graphic  location  for  the  parameters  and  categories 
which  they  tested,  "Why  should  low  ceilings  constitute 
an  exception?"  The  answer  appears  to  lie  in  the  non¬ 
homogeneity  of  conditions  which  can  produce  low  ceilings 
compared  to  those  less  readily  affected  by  the  radia- 
tional  transfer  process  occurring  near  the  atmospheric 
lower  boundary.  Thus,  the  correlation  may,  in  fact,  be 
quite  constant  for  each  type  of  ceiling-producing  cloud. 
Since  fog  and  stratus  (the  main  cause  of  low  ceilings) 
are  highly  dependent  on  the  local  solar  hour,  steps  were 
taken  to  partition  the  data  according  to  the  local  time 
of  day. 
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To  determine  the  optimum  partitioning,  we  decided 
to  assume  that  a  generalized  solution  to  the  problem 
of  determining  conditional  probabilities  statistically 
(once  their  unconditional  counterparts  were  derived 
observationally)  does,  in  fact*  exist  providing  the 
hourly  and  monthly  variability  of  cloud  types  are  con¬ 
strained.  An  "inexhaustible"  data  base  was  devised  for 
testing  the  array  of  options  open  to  us  for  finding 
that  solution  by  arbitrarily  selecting  and  combining  the 
hourly  station  records  of  fifteen  AWS  widely  separated 
locations  (see  Figure  5).  The  respective  smoothed 


unconditional  frequency  of  occurrence  and  their  2-hour 
recurrence  probabilities  for  every  hour-of-the-day  and 
month-of-the  hear  were  extracted  for  each  particular 
station.  A  common  diagram  was  then  produced  using  this 
set  of  statistics  for  each  of  the  stations  assuming 
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each  hour-of-the-day  constituted  an  initial  forecast 
hour  (tc).  The  unconditional  frequencies  of  low 
ceiling  occurrences  at  forecast  time,  tQ»  were  repre¬ 
sented  on  the  ordinate.  The  corresponding  unconditional 
values  2  hours  later  define  the  abscissa.  The  respec¬ 
tive  2-hour  recurrence  probabilities  were  plotted  at 
the  coordinate  intersection  of  these  two  values.  Since 
all  these  informations  were  derived  from  observed  values, 
their  plotted  array  represents  a  host  of  observed  ceil¬ 
ing  occurrences  and  recurrences  under  widely  variant 
topographic  and  geographic  conditions.  As  expected, 
when  all  possible  times  of  the  day  and  months  of  the  year 
were  included  in  the  plotted  array,  an  analysis  of  the 
recurrence  statistics  was  difficult  to  complete  under 
required  tolerances.  However,  as  more-and-more  con¬ 
straints  were  imposed,  such  that  only  those  select  por¬ 
tions  of  each  station's  data  base  which  pertains  to 
certain  solar  times  and  seasons  were  plotted  on  the  common 
graph,  the  variance  in  the  recurrence  data  about  the 
analyzed  value  became  progressively  reduced.  Essentially 
these  results  confirm  our  original  assumption  that  when 
the  ceiling-producing  clouds  are  rendered  near  homogene¬ 
ous  through  considerations  based  upon  the  time  of  day  and 
season  of  the  year,  the  concepts  of  universality  as  ad¬ 
vanced  by  Gringorten  essentially  pertain  for  the  low 
ceiling  category  as  well.  Hence,  if  two  or  more  stations 
exhibit  identical  profiles  of  unconditional  probabilities 
within  any  segment  of  the  time  partitioned  subset,  the 
extent  to  which  they  were  topographically  or  geographically 
induced  becomes  essentially  irrelevant  for  our  purposes. 
Thus,  a  station's  records  contains  at  least  two  types  of 
pertinent  information  to  the  problem  at  hand.  Locally 
dependent  geographical  and  topographical  factors  work  to 
identify  which  respective  coordinate  intersections  apply 
to  each  particular  station.  The  recurrence  statistic  at 
each  given  intersection  does  not  require  a  knowledge  of 
which  stations  meet  the  particular  magnitude-and- slope 
criteria  for  their  frequency  of  occurrence  values  which 
are  demanded  by  that  specific  coordinate  intersection. 

After  considerable  experimentation,  it  was  found  that 
a  16-fold  partitioning  of  the  data  (four  6-hour  periods 
per  day  for  each  season  of  the  year)  defined  relation¬ 
ships  between  event  occurrence  frequencies  and  their  re¬ 
currence  probabilities  which  were  sufficiently  independent 
of  geography  and  topography  that  these  features  need  not 
be  identified  more  specifically.  Two  such  graphs  are 
shown  in  Figures  6a  and  6b.  The  diurnal  difference  in 
these  relationships  is  seen  by  comparing  the  isopleth 
slopes  between  the  two  graphs.  Complete  sets  of  2-  and 
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fig.  I«.  7-hour  occurrence- recurrence  relet ionships 
deduced  fro*  the  l$-3tetien  dete  bene 
during  the  evening  hourt  (1700-23Q0L)  for 
Sueuncr  ( June- August ) . 


Fig.  tb.  7-hour  occurrence-recurrence  relationships 
deduced  fro«  the  lS-steticn  det*  ben* 
during,  the  night  hours  <  ?3‘)0-tS00L)  for 
Suweer  (June-August). 
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4-hour  ceilings  and  2-hour  visibility  graphs  are  avail¬ 
able  in  a  Saint  Louis  University  publication,  Hull  (1973)* 
They  represent  another  by-product  of  our  continuing  re¬ 
search  efforts.  Essentially,  they  provide  a  means  of 
generating  conditional  probabilities  for  those  stations 
having  weather  records  of  too  short  a  time  duration  to 
permit  their  determination  by  direct  computer  tracing. 

The  accuracy  of  the  graphical  relationships  was 
tested  against  dependent  and  independent  data.  The  de- 
pendent-data  verifications  were  accomplished  by  simply 
using  the  graphs  and  the  frequency-of -occurrence  data 
to  reconstitute  the  2-hour  recurrence  arrays  for  each 
of  the  15  stations  in  formats  such  as  that  shown  in 
Figure  3.  These  were  then  compared  against  the  AWS 
generated  and  objectively  smoothed  recurrence  proba¬ 
bilities  for  each  individual  station,  category  and  time 
period.  These  comparisons  readily  verify  the  fact  that 
the  techniques  developed  did,  in  fact,  apply  to  each 
individual  station  as  well  as  to  the  majority.  Figures 
7  and  8  show  the  graphically  reproduced  arrays  for 
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Fig.  7.  Estimated  2-hour  recurrence  probabilities 

for  ceilings  equal  to  or  less  than  200  feet 
for  Randolph  AFB ,  Texas. 
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Fig.  8.  Estimated  2-hour  recurrence  probabilities  for 
ceilings  equal  to  or  less  than  200  feet  for 
Vandenberg  AFB,  California. 

two  of  the  fifteen  stations.  Their  observed  counterparts 
are  shown  in  Figures  3  and  9,  respectively.  These  two 

*0tw  *  *  *  «  ■  in  i?  i*  m  »•  ?n  ?? 


Fig.  9.  Observed  conditional  2-hour  recurrence  proba¬ 
bilities  for  ceilings  equal  to  or  less  than 
200  feet  for  Vandenberg  AFB,  California. 
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THE  J-HR  RECONSTRUCTED  CONDITIONAL  PROBABILITY  Of  CriLINGS«J00'  RCCURRINO 


Elf.  11.  March  ATS,  Calif. 


Again, the  estimated  2-hour  recurrence  probabilities  were 
obtained  by  introducing  observed  frequency-of-occurrence 
pairs  for  each  of  the  three  respective  stations  into  the 
appropriate  graphs.  The  recurrence  statistics,  provided 
by  the  AWS  computer-tracing  efforts,  against  which  these 
graphical  estimates  are  to  be  compared,  are  shown  in 
Figures  13,  14  and  15. 


THE  2-HR  OBSERVED  CONDITIONAL  PROBABILITY  OF  CEILINGS<  200'  RECURRING 
HOUR  0  2  H  6  I  10  12  14  16  II  20  22 


70 


Fig.  13.  Wiesbaden,  Germany 


18 


fif.  15.  March  AFB.  Calif 
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The  same  general  methodology  as  that  presented 
above  for  the  0-2  hour  time  period  also  was  found  to 
pertain  for  the  4-hour  time  frame.  Likewise,  it 
proved  valid  when  visibilities,  rather  than  ceilings, 
represented  the  parameters  being  investigated.  These 
graphical  relationships  are  similarly  contained  in 
the  aforementioned  Saint  Louis  University  publication. 

III.  A  REVISED  MODEL  WHICH  INCORPORATES  THE  LATEST  OBSERVED 
WIND  INFORMATION 

The  fact  that  the  persistency  of  a  ceiling  event 
is  dependent  on  the  wind  regime  occurring  in  associ¬ 
ation  with  it  is  shown  in  Figure  16  where  the  recur¬ 
rence  data  for  Randolph  AFB,  Texas  are  presented  as  a 
mere  example . 


0  2  *  <>  *  10  1*  14  16  J®  2n  a2 


Fig.  16.  Estimated  probabilities  that  observed  ceilings 
equal  to  or  less  than  200  feet  .  presently 
occurring  within  a  given  wind  sector,  will  re¬ 
occur  2  hours  later  for  Randolph  AFB  during 
the  period  December  through  March. 


It  is  readily  evident  from  these  data  that  the  problem 
is  not  one  of  whether  or  not  wind  information  should  be 
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incorporated  into  the  persistency  evaluations  but 
rather  one  of  how  to  effect  it  without  compromising 
the  general-applicability  attributes  of  the  modelled 
products.  Two  other  pertinent  questions  also  arise. 

How  severely  can  a  given  station's  data  base  be  par¬ 
titioned  and  still  be  expected  to  produce  climatically 
significant  conditional  probabilities?  Secondly,  do 
the  empirical  relationships  which  we  have  found  still 
pertain  when  stratified  subsets  of  a  station's  data 
base  are  considered?  Neither  of  these  questions 
afford  an  obvious  answer.  The  number  of  years  of  his¬ 
torical  data  required  to  establish  climatic  signifi¬ 
cance  for  the  computed  results  depends  upon  the  per¬ 
centage  frequency  that  an  event  occurs  and  the  variance 
in  its  persistency  qualities.  Often  both  of  these  are 
enhanced  by  dividing  the  data  into  subsets.  Thus,  the 
question  of  the  number  of  years  needed  to  establish 
climatic  significance  for  a  given  type  is  not  a  simple 
multiple  of  that  needed  for  the  unstratified  data  base. 
This  point  becomes  particularly  evident  when  moisture 
considerations  are  invoked  as  the  stratifying  agent. 

Since  certain  humidity  values  preclude  the  formation 
of  clouds,  and  others  almost  guarantee  it,  this  type  of 
stratification  drastically  alters  the  frequency  of 
occurrence  statistics.  With  respect  to  the  question  as 
to  whether  or  not  the  empirical  relations  given  by  the 
graphs  apply  to  stratified  subsets,  the  answer  is  "yes" 
providing  the  stratifying  agent  itself  persists.  This 
conclusion  was  determined  by  introducing  the  3 -wind 
sector  stratified  "unconditional"  probabilities  con¬ 
tained  in  the  AWS  statistics  into  the  graphs  of  Section  II 
and  comparing  the  estimated  recurrence  probability  values 
against  their  computer-traced  counterparts. 

Certain  problems  also  arise  in  using  these  graphs 
on  the  stratified  subsets  which  are  not  found  when  un¬ 
stratified  climatologies  are  considered.  The  problem  is 
simply  one  of  determining  which  wind  sector  to  use  for 
obtaining  the  wind-stratified  frequency  of  occurrence 
statistic  for  the  hour  tQ  +  2.  This  obviously  requires 
a  2-hour  forecast  of  the  winds  and  an  assumption  that  a 
changing  wind  direction  does  not  represent  a  sufficiently 
significant  change  in  "physics"  to  destroy  the  Markov- 
chain-type  relationships  in  the  data.  Since  low  ceilings 
frequently  occur  or  persist  in  association  with  weak 
winds  and  dissipate  when  the  winds  accelerate,  particu¬ 
larly  in  the  early  morning  hours,  the  graphs  were  found 
to  be  biased  towards  persistency.  Herein  lies  the  most 
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formidable  obstacle  we  have  faced  in  devising  a  sta¬ 
tistical  model  which  brings  additional  meteorological 
knowledge  to  bear  upon  the  persistency  of  an  event. 

The  approach  we  are  using  demands  that  the  same  know¬ 
ledge  be  available  for  tQ  +  2  as  that  used  to  stratify 
the  data  at  tQ.  Thus,  care  must  be  taken  that  one 
forecast  problem  is  not  simply  being  substituted  for 
another.  The  statistics  for  Randolph  and  Scott  AFB 
shown  in  Figure  17  emphasize  this  problem. 


3  Wind  3  Wind 

Sector  Sector 

Strati-  Strati- 
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Persist¬ 
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56 
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69 

Fig.  17.  2-hour  forecast  verification  for  cases  when  the 
ceiling  initially  ranged  between  0  and  199  feet 
(winter  season). 


Verification  statistics  for  ceilings  and  winds  were  ex¬ 
tracted  from  the  observed  records  for  each  station  when¬ 
ever  the  time-marched  initial  time,  to,  showed  ceilings 
less  than  200  feet  to  be  occurring.  The  appropriate 
3  wind-stratified  frequency  of  occurrence  statistics  for 
determining  the  ordinate  position  on  the  graphs  were  then 
extracted  from  presentations  such  as  those  given  in 
Figure  1.  Two  sets  of  abscissa  values  were  also  selected. 
One  neglected  wind  information  throughout  the  two  hour 
time  period.  The  other  was  based  on  a  knowledge  of  the 
wind  initially  and  two  hours  later.  Whenever  the  esti¬ 
mated  conditional  probabilities  for  tQ  +  2  were  50%  or 
greater,  a  persistent  forecast  was  rendered.  The 
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persistency  of  low  ceilings  during  the  3  year  winter  period 
for  both  stations  is  greater  than  50%  for  the  2  hour 
period  (Figure  17,  Col.  2).  Columns  3  and  4  show  the 
percentage  of  correct  forecasts  under  all  wind  conditions, 
while  columns  5  and  6  show  the  results  when  wind  strati¬ 
fications  were  invoked.  Note  that  the  graphs  of 
Section  II  faithfully  provided  verification  statistics  com¬ 
parable  to  those  given  by  the  AWS  computer-traced  product 
for  both  locations  when  wind  information  was  not  con¬ 
sidered,  i.e.,  when  the  model  discussed  in  Section  II  was 
used.  Also  note  that  the  estimated  conditional  proba¬ 
bilities  verified  better  for  the  wind-stratified  data  for 
Randolph  than  they  did  when  wind  information  was  not  con¬ 
sidered.  The  data  for  Scott  showed  different  results. 

Here  wind  stratification  proved  less  effective  due  to  the 
simple  fact  that  low  ceilings  at  this  mid-continent  lo¬ 
cation  are  more  critically  affected  by  a  more  variable 
wind  climatology. 

Next  the  emphasis  was  shifted  from  one  of  low  ceil¬ 
ings  persisting  to  that  of  low  ceilings  forming  rapidly 
from  a  much  improved  height  category,  i.e.,  the  rapidly 
deteriorating  weather.  The  records  of  Randolph  AFB  were 
carefully  screened  to  separate  out  those  cases  when  rapid 
ceiling  deteriorations  occurred  and  the  data  were  proc¬ 
essed  through  the  appropriate  graphs.  Some  case  studies 
are  shown  in  Figures  18  and  19. 
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Fig.  1*. 


"Unconditional"  probability  of  avant 
occurrence*  at  Randolph  AFB  on 
December  23,  19S1. 
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(whan  neither  latest  observed  dew-point  depressions  or  winds  were 
considered) . 
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Fig.  1*.  "Uncondit ional"  probability  of  event 
occurrences  *t  Randolph  AFB  on 
January  IS ,  1919. 


Note  that  neither  the  wind  stratified  climatic 
statistic,  nor  those  devoid  of  such  considerations  pro¬ 
duced  conditional  probabilities  for  any  category  other 
than  "F"  remotely  approaching  the  50%  criterion  which 
we  arbitrarily  set  for  determining  that  a  new  forecast 
be  rendered.  Hence,  climatic  considerations  based  on 
wind  fields  alone  appear  to  offer  little  encouragement 
to  the  problem  of  forecasting  rapid  deteriorations. 

The  limited  data  does,  however,  suggest  that  considera¬ 
tions  based  upon  moisture  fields  (note  that  the  tempera¬ 
ture  dew-point  spread  indicates  increasing  humidities 
long  before  the  ceiling  is  found  to  lower  or  form)  may 
deserve  a  better  fate.  This  point  was  independently 
supported  by  each  of  the  AWS  teams  which  attempted  to 
use  our  products  on  a  trial  basis.  It  represents 
another  topic  which  will  be  presented  later  in  the  report. 
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Before  proceeding  to  that  discussion,  however,  it  should 
be  noted  that  the  wind-stratified  statistics  provide 
much  useful  information  apart  from  that  of  statistically 
deducing  recurrence  probabilities. 

Presentations  such  as  that  shown  in  Figure  20 
quickly  alert  the  forecaster  to  those  ceiling  forecast 
problems  at  his  location  which  are  the  most  strongly 
wind  related.  Of  particular  significance,  these  statis¬ 
tics  are  not  predicated  on  a  knowledge  of  the  observed 
ceiling  conditions  at  forecast  time.  They,  therefore, 
apply  equally  well  to  tomorrow's  weather,  for  example, 
as  they  do  for  today's.  Hence,  they  can  be  used  in  con¬ 
junction  with  the  numerical  models  for  any  time  period 
by  simply  asking  that  model  to  supply  a  gross  estimate 
of  the  wind  direction  when  the  speeds  exceed  a  certain 
limit.  We  used  4  mph  as  that  criterion  in  our  own  study. 
For  example,  if  the  24-hour  fine  mesh  upper  air  and  sur¬ 
face  progs  show  a  low-level  wind  out  of  the  east  at  0600 
LST,  the  forecaster's  odds  that  ceilings  £199  feet  will 
form  are  given  by  Figure  20  to  be  0.157.  The  probabilitv 


Local  time  of  day 


Fig.  20.  Observed  frequencies  for  ceilings  equal  to  or 

le68  than  200  feet,  occurring  within  given  wind 
sectors,  for  Randolph  AFB  during  the  period 
December  through  March. 
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two  hours  later,  providing  the  wind  is  forecast  to 
persist,  is  found  to  be  0.120.  The  2-hour  persis¬ 
tency  probabilities  for  this  forecast  ceiling  con¬ 
dition  starting  at  0600  can  be  computed  from  the 
graphs  to  be  about  0.59  (59%)(see  Figure  16. 

Such  applications  of  the  wind-stratified  tech¬ 
nique  should  only  be  attempted  when  the  wind  direction 
can  be  assumed  to  persist.  The  reason  being  that  a 
shifting  wind  often  infers  a  discontinuity  in  the 
physical  processes  that  maintain  low  ceilings.  For 
this  reason  we  turned  our  attention  to  additional 
stratifying  agents  of  a  more  continuous  nature  than 
that  of  winds  alone.  These  will  be  described  in  the 
section  to  follow. 


IV.  A  REVISION  TO  INCLUDE  MOISTURE  CONSIDERATIONS  (TEMPERATURE 
DEW-POINT  SPREADS) 

Some  problems  arising  in  conjunction  with  the  model¬ 
ing  of  wind-stratified  data  to  provide  recurrence  proba¬ 
bilities  were  discussed  above.  Others  arise  from  the 
fact  that  winds  are  not  physically  related  to  ceilings 
but  essentially  enter  the  picture  by  providing  the 
mechanisms  which  amass  those  parameters  that  are  critical 
to  ceiling  informations.  Moisture,  on  the  other  hand, 
plays  a  more  direct  role  since  it  is  a  fundamental  in¬ 
gredient  in  cloud-droplet  formation  processes. 

For  this  reason  humidity  parameters  (we  used  the 
temperature  dew-point  spread)  appear  to  be  superior  to 
winds  as  a  stratifying  agent  since  an  imposed  humidity 
constraint  more  effectively  serves  to  reduce  the  variances 
in  the  occurrence  and  recurrence  statistics.  This  offers 
the  very  fundamental  advantage  that  climatically  signifi¬ 
cant  statistics  can  be  realized  from  a  smaller  period  of 
records.  Another  advantage  is  that  certain  problems 
arising  with  the  rarely  occurring  ceiling  event  can  be 
effectively  recast  to  that  of  a  rarely  occurring  continuous 
moisture  parameter.  This  also  has  its  benefits  by  per¬ 
mitting  a  maximum  of  information  to  be  extracted  from  a 
minimum  of  data.  To  illustrate  these  points,  an  inspec¬ 
tion  of  the  historical  records  readily  shows  that  ceilings 
^199  feet  occur  rather  infrequently  for  most  stations. 
Hence,  the  statistics  are  often  computed  on  too  small 
a  data  sample  for  them  to  stabilize  due,  in  part,  to  the 
variability  of  the  recurrence  statistics.  Note,  for  ex¬ 
ample,  in  Figure  1  during  the  daylight  hours  of  the  winter 
season  unconditional  frequencies  are  found  to  be  as  low  as 
.008.  Figure  21  shows  this  same  statistic  to  be  much 
higher  when  ceiling  occurrences  are  "conditioned"  on  the 
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Fig.  21.  Probability  (percent)  of  Cloud  A  as  a  function 
of  dew  point  (all  winds)  at  Randolph  AFB,  Texas 
in  the  winter  season. 


associated  temperature  dew-point  spread.  The  argument 
is  not  one  that  by  some  wizardry  a  rare  occurring  event 
now  becomes  a  highly  occurrent  one.  Rather,  the  criter¬ 
ion  for  determining  the  rarity  is  shifted  to  the  occur¬ 
rence  of  a  temperature  dew-point  spread.  Although  these 
spreads  rarely  occur,  the  statistics  associated  with 
them  can  better  be  inferred  by  bringing  the  entire  data 
records  to  bear  upon  the  problem  as  discussed  in  con¬ 
junction  with  Figure  1.  The  moisture  parameter  affords 
other  distinct  advantages  which  soon  become  clearly  evi¬ 
dent.  For  example,  each  specific  ceiling  category  can 
be  represented  by  some  median  height  value  within  that 
category.  At  this  cloud  base  height  the  humidities  are 
automatically  pegged  as  being  near  100%,  otherwise  a 
ceiling  would  not  exist.  This  knowledge  of  height  and 
humidity,  together  with  that  of  moisture  values  at  the 
surface,  establishes  a  linear  approximation  to  the 
moisture  lapse  rate  between  the  ground  and  the  ceiling 
base.  This  information  provides  a  crude  "typing"  of 
the  ceiling-producing  clouds.  Primarily  it  serves  to 
delineate  the  low-level  fog  cases  from  the  remainder  of 
the  sample.  The  fact  that  fog  has  different  persistency 
characteristics  than  other  types  of  ceiling-producing 
phenomena  is  common  knowledge . 

Displays  such  as  those  shown  in  Figure  22  were  pro¬ 
duced  for  each  ceiling  category.  Note  how  the  frequency 
of  ceiling  occurrences  vary  diurnally  with  surface 
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Fig.  22.  Ceiling  frequencies  of  clouds  1,000  to<3,000  feet  as  a 
function  of  dew-point  depression  and  time  of  day  for 
Randolph  AFB  (winter). 


humidities.  Although  similar  variations  are  evident  for 
other  ceiling  categories  as  well,  this  particular  illus¬ 
tration  was  chosen  since  these  ceilings  lie  in  the  upper 
reaches  of  the  atmosphere's  lower  boundary  layer  and  have 
the  option  of  rising  above  it  or  deteriorating  more  deeply 
into  it.  The  diurnal  relationships  between  ceilings  and. 
given  temperature  dew-point  spreads  presumably  reflect  diur¬ 
nal  changes  in  low-level  thermal  stratifications  associated 
with  the  local  solar  hour.  Thus  another  bit  of  pertinent 
information,  having  a  bearing  on  ceiling  persistency 
appears  to  be  implicitly  incorporated  into  the  model. 
Moisture  stratified  subsets  of  the  data  for  a  number  of 
stations  were  processed  to  produce  similar  relationships 
to  those  shown  in  Fig*  6a  and  6b  between  moisture  stratified 
occurrence  values  and  the  recurrence  probabilities  that 
a  ceiling  occurring  under  these  conditions  will  recur  two 
hours  later.  Matrices  of  information  were  extracted  from 
each  of  these  respective  graphs  and  placed  on  punch  cards 
and  magnetic  tape  to  facilitate  the  calculation  of  re¬ 
currence  probability  statistics  for  any  given  station  and 
ceiling  event  (see  Figure  23). 


28- 


|*IHa|.  »*t>  >0|f,T  b»»r*V  0  IN|T1,1  CL«‘H>  A  ST*IfO«  «lNo  CALM  Al*TlP 

!•  z  HmiU  C0»4^rT|0HAL  *»POHAd|LlT¥ 

PINAL  riH*L  *'  r  LOCAL  ^T«Nf)«**0  ll-t 


JPKtAO 

CLULp 

1 

3 

A 

.  S 

6 

7 

4 

9 

10 

U 

12 

13 

1* 

IS 

1 6 

It 

id 

19 

2o 

2l 

22 

23 

24 

..  A 

•  9  . 

«4 

99 

92 

•  o 

•1 

do 

7« 

71 

72 

77 

77  . 

76 

77 

70 

•  *- 

.  to 

7d 

7A 

74 

75 

79 

93 

-*3 

9 

7 

5 

6 

9 

11 

11 

12 

1* 

2> 

21 

19 

13 

44 

16 

IS 

1*1 

-l? 

17 

13 

15 

45 

11 

A 

A 

c_ 

7  > 

z 

,  2 

2  _ 

6 

.4 

4. 

_  s 

_ 1 

X— 

4— 

4 - 

4- 

1 

—  3  • 

- *- 

U 

1 

1 

2 

1 

2 

1 

1 

1 

1 

1 

1 

1  ~ 

1 

1 

i 

2 

2 

’  2 

1 

1 

1 

1 

1 

..  t 

1 

1 

1 

1 

1 

1 

1 

1 

.1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

..  I  . 

1 

•  1 

1 

1 

f 

S 

5 

A 

5 

A 

A 

2 

1 

0 

0 

0 

0 

0 

0 

0 

1 

2 

3 

5 

1 

A 

4 

4 

4 

A 

s* 

5l 

52 

A9 

ah 

41 

51 

50 

*d 

46 

A4 

64 

9? 

♦  3 

95 

46 

44 

45 

44 

A6 

AO 

Si 

55 

55 

b 

21 

2? 

1" 

19 

19 

29 

29 

3*_ 

90_45 

*7_ 

.45  _ 

*7 

AS 

9ft 

3« 

3a 

31 

27_ 

it.. 

25 

25 

23. 

— XL— 

1 

C 

10 

In 

if 

19 

19 

13 

12 

9 

6 

7 

6 

6 

d 

H 

19 

h 

in 

19 

13 

12 

11 

IP 

9 

11 

u 

s 

9 

3 

2 

S 

A 

2 

1 

1 

1 

1 

1  . 

2 

3 

2 

3 

3 

6 

6 

.  4 

A 

S 

4 

...J.. 

t 

1 

1 

2 

A 

1 

1 

1 

1 

2 

1 

n 

1 

1 

1 

2 

l 

> 

1 

1 

2 

2 

1 

1 

1 

F 

9 

In 

9 

9 

9 

7 

5 

3 

1 

0 

.  0 

0. 

0 

0 

1 

1 

S 

9  9  . 

-  4 

• 

d  . 

...  ft  . 

A 

2r 

19 

17  . 

17 

17 

1  9 

19 

24  _ 

19 

14 

1 6  _ 

.  17_. 

.16_ 

16 

16 

15-45 

15 

16-.I7-. 

24—  23 

23—21— 

• 

*0 

if) 

2? 

2n 

2 1 

?  7 

72 

*7 

6  0 

66 

64 

69 

61 

Ss 

5ft 

A  1 

31 

?• 

25 

24 

22 

2n 

Id 

20 

i 

c 

39 

37 

3S 

19 

1« 

23 

23 

?* 

19 

11 

In 

10. 

16 

23 

25 

29 

34 

3* 

34 

33 

34 

3a 

37 

37.- 

u 

5 

9 

A 

23 

2A 

1H 

17 

J 

3 

2 

f 

3 

3 

i 

S 

6 

Q 

6 

A 

7 

A 

6 

5 

5 

L 

1 

1 

1 

2 

3 

s 

2 

2 

2 

2 

\ 

1 

1 

2 

2 

2 

1 

.  5 

1 

-  1^ 

1 

1 

1 

-  1 

» 

1* 

19 

19 

in 

1  7 

9 

7 

* 

2 

1 

n 

0 

1 

1 

2 

S 

T 

1« 

id 

Id 

17 

16 

1ft 

1* 

A 

9~  " 

M 

% 

A 

9 

*9 

”,0 

*d” 

’  6 

4 

A 

3 

~3~ 

~3~ 

1 

4 

4 

4 

5 

S 

7 

9 

'• 

e 

17 

17 

16 

lA 

IS 

2« 

2* 

M 

Afo 

ss 

6, 

55 

92 

3d 

d| 

91 

19 

1* 

15 

13 

lA 

13 

11 

15 

3 

C 

19 

19 

15 

IS 

1* 

1» 

30 

13 

26 

24 

2a 

29 

Aft 

*2 

9ft 

Si 

29 

t* 

16 

19 

if 

26 

?• 

22 

L 

32 

3 1 

31 

11 

12 

11 

23 

1* 

IS 

11 

7 

in 

12 

13 

9 

17 

29 

35 

3d 

3b 

33 

29 

2d 

29 

L 

1 

1 

1 

19 

21 

?3 

S 

2 

2 

2 

? 

l 

2 

2 

A 

A 

• 

1 

1 

1 

1 

1 

1 

1 

- 

_  .  f 

29. 

27 

29 

33 

-3-j 

19 

9 

• _  J 

1 

- 1- 

--1- 

1 

-  2 

3 

7 

•  11 

2* 

24. 

26 

25 

24 

23 

-25- 

A 

3 

A 

A.. 

~.J 

.  3 

3. 

.  3 

A 

.  3 

.2. 

— -1- 

..  A.. 

0 

..  0  .. 

0 

i 

...a- 

~1~ 

-.X.. 

-1.. 

...  2 . 

b 

? 

9 

in 

9 

9 

11 

13 

1 5 

16 

2  A 

?A 

27 

2n 

15 

1? 

1  <* 

H 

4 

7 

7 

6 

5 

4 

4 

A 

C 

15 

If 

15 

IS 

1A 

17 

26 

A* 

52 

53 

SS 

ST 

59 

*S 

?S 

n 

id 

12 

12 

15 

16 

17 

u 

15 

19 

12 

13 

15 

15 

14 

2* 

19 

16 

ll 

15 

19 

23 

J* 

A  7 

5> 

A/ 

AA 

A3 

93 

43 

17 

1  7 

_ 

L 

IV 

19 

19  _ 

2n 

2a 

21 

22 

,S._ 

-X. 

3  - 

- \. 

- 1— 

?  * 

-2- 

4— 

— 1 — 1—  1- 

1 

-  2ft - 

F 

♦  1 

*0 

*1 

An 

3v 

33 

20 

/ 

s 

2 

? 

2 

2 

3 

6 

9 

26 

31 

35 

36 

3a 

33 

3ft 

40 

A 

1 

1 

Y‘ 

i 

1 

1 

1 

1 

1 

0 

n 

0 

0 

U 

0 

0 

0 

0 

0 

0 

•9 

0 

1 

1 

• 

• 

5 

7 

6 

6 

7 

M 

9 

9 

in 

7 

ft 

6 

S 

A 

3 

3 

i 

2 

2 

2 

1 

1 

2 

s 

c 

11 

1* 

9  * 

11 

1? 

19 

19 

3« 

-.1 

47 

66 

63 

S* 

36 

19 

11 

4 

s 

6 

7 

7 

4 

11 

11 

0 

19 

17 

19 

1A 

17 

15 

17 

3d 

26 

27 

24 

.33 

S  V  . 

62 

6ft 

_ S  7 

51 

46 

_  aS 

♦  S 

25 

22 

20 

c 

23 

2? 

21 

13 

IS 

14 

Id 

2 

A 

2 

2 

3 

3 

7 

a 

1 

3 

2* 

A 

3 

23 

22 

22 

F 

*S 

*5 

99 

55 

Sa 

47 

37 

|3. 

7. 

3 

.  .3 

... X- 

...  A.. 

« 

.4*. 

3  1 

3* 

2X 

.42 

*3 

♦  ?. 

43 

..  A 4 

A 

9 

.  1 

1 

1 

1  . 

1 

1 

0 

0 

...  0. 

...0 

0 

.  0 

0 

0 

...ft 

.  ..2 

.  0.. 

_ ft 

.  0- 

...  0  . 

-.4. 

? 

A 

A 

3 

A 

A 

6 

6 

6 

3 

2 

2 

? 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

_  6 

_ C 

4 

9 

..  v 

9 

H 

9 

9 

_  3o 

36 _ 2, 

1 6._, 

.1,1 

.  A 

J- 

i _ 1 

_ X  _ 

_ A^_ 

_  A 

_  fc_— 2. 

"  1»  19 

15 

”l  9 

19  15 

7«' 

14 

37 

A  1 

5? 

SS 

65 

b9 

9) 

3’ 

36 

i* 

25 

23 

21 

If 

2* 

20 

c 

25 

2a 

-  13 

12 

1? 

13 

id 

12 

23 

IS 

2 

2 

.  3 

7 

49 

24 

25 

29 

27 

16 

15 

14 

to 

25 

f 

♦  7 

♦  7 

59 

M 

41 

44 

Ad 

3* 

9 

6 

5 

S 

6 

7 

16 

27 

34 

a! 

44 

41 

59 

5* 

59 

4ft 

H 

2 

-  3 

.  4. 

3 

*3 

3 

*  . 

* 

_ 2 

_ 1 

-  L 

.  1 

1. 

~  l 

^i- 

_ 1 

—  1 

„|. 

— 1 

-  1- 

.  -4. 

“  r  “ 

C 

7 

7 

6 

6 

6 

7 

6 

'“*V 

\* 

16 

4 

6 

m 

1 

7 

1 

1 

l 

2 

3 

1 

4 

A 

V 

1* 

IS 

13 

13 

13 

13 

3* 

«d 

M 

67 

«2 

66 

51 

*3 

36 

J9 

2' 

ll 

id 

U 

IT 

)f 

Id 

t 

11 

U 

17 

11 

11 

1* 

IS 

ftft 

21 

u 

l 

21 

34 

JT 

Jft 

*9 

IV 

17 

15 

14 

17 

10 

U 

__  _ f 

_44  . 

9* 

*5 

•  7 

•  7 

63 

56 

3* 

*3 

.  T. 

..A. 

,..6 

X. 

..  d 

16 

31- 

.30 

S± 

AO 

.44 

.66  .61 

46 

-.44. 

9 

1 

2 

3 

3 

2 

2 

2 

4 

1 

1 

1 

J 

1 

1 

1 

1 

1 

l 

1 

\ 

1 

1 

1 

« 

C 

s 

S 

5 

A 

S 

6 

4 

a 

4 

4 

1 

3 

2 

1 

1 

.1 

.  1 

1 

I 

.  1 

2 

3 

A 

0 

19 

1A 

12 

11 

1 0 

n’ 

1« 

19 

*1 

S* 

7  1 

72 

At 

Su 

'2 

36 

2* 

22 

Li 

1* 

12 

14 

19 

Id 

t 

9 

In 

10 

l  n 

11 

i? 

13 

1* 

22 

17 

14 

11 

26 

)* 

JA 

In 

21 

id 

13 

1* 

15 

12 

11 

10 

r 

72 

72 

69 

41 

*s 

3  a 

19 

7 

6 

T 

9 

22 

it 

A  T 

47 

67 

79 

'0 

70 

4* 

6d 

*  ■ 

A' 

~r 

~’o 

— 

~  0 

0 

— i 

«  1 

0- 

0 

0 

" 

0 

"'fl 

a 

"V 

0 

«“ 

*" 

i 

4 

0 

D 

4 

ft 

i 

2 

2 

2 

2 

2 

t 

1 

1 

1 

1 

1 

1 

1 

l 

1 

l 

i 

1 

x 

1 

j 

1 

1 

9 

c 

9 

S 

5 

A 

A 

S 

5 

s 

A 

2 

7 

l 

i 

1 

l 

1 

1 

i 

I 

k 

? 

J 

3 

4 

U 

12 

1? 

11 

9 

9 

9 

15 

J9 

A» 

64 

66 

6? 

Ad 

Jd 

31 

?1 

1 7 

32 

6 

11 

11 

IS 

l* 

t 

9 

9 

10 

12 

12 

13 

1* 

21 

?* 

t> 

25 

29 

40 

Jt 

3” 

2* 

i» 

1* 

IS 

11 

If 

ft 

10 

— 

f 

_ 7  2 

_7a_ 

J!  5 

7a 

72 

*5 

5*»__35 

25 

16 

7 

LU 

26 

-31 

ll 

.  A* 

72_ 

_n 

75 

73 

72 

-71- 

¥ 

1 

1 

2 

1 

1 

1 

1 

1 

1 

l 

1 

1 

l 

1 

1 

1 

1 

I 

1 

i 

1 

1 

l 

1 

- 

C 

3 

3 

3 

3 

3 

A 

3 

J 

2 

I 

1 

l 

1 

1 

i 

1 

1 

3 

1 

i 

l 

r 

3 

A 

D 

13 

13 

9 

7 

7 

1  n 

IS 

16 

*1 

4  6 

53 

5n 

A* 

ii 

?4 

17 

ll 

T 

ii 

9 

17 

13 

H 

L 

9 

A 

9 

11 

12 

12 

1* 

|  9 

22 

25 

3.1 

40 

40 

Ja 

2* 

19 

1* 

1 3 

* 

13 

1  * 

9 

9 

f 

7  A _ 7^ 

77 

79 

7  7 

73 

hi 

4i 

39 

12 

in 

t 

4 

ll 

4‘* 

A  4 

At 

71 

76 

if 

76 

7% 

74 

7* 

Fig,  23  Computer  produced  2-hour  conditional  probabilities 
using  dew-point  spread  and  wind  stratified 
"universal"  relationships  for  Randolph  AFB,  initial 
spread  0,  for  initial  ceilings  less  than  200  feet, 
(Note:  the  letters  "A , B ,C ,D  ,E ,F"  refer  to  cloud 
categories  used  by  AWS*  The  definitions,  in  order, 
are  "ceilings  0  to  <200  feet,  200  to  <500  feet, 

500  to  <1000  feet,  1000  to  <3000  feet,  3000  to 
<10000  feet,  and  ceilings  at  or  above  10000  feet  or 
no  ceiling.") 
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A  major  obstacle  common  to  that  found  in  Section  III 
remains  with  us,  however.  How  is  one  to  determine  the 
dew-point  spread  two  hours  hence?  This  is  a  fundamental 
requirement  of  the  technique.  Two  possible  solutions 
came  to  mind.  One  was  to  provide  the  forecaster  a  clima¬ 
tology  of  dew-point  spread  changes.  This  climatology  was 
obtained  by  searching  the  historical  records  and  noting 
how  similar  initial  moisture  conditions  subsequently 
evolved  in  the  next  two  hours.  This  approach  was  quickly 
shown  to  be  limited  (both  in  our  evaluations  and  those  of 
the  field  stations)  simply  because  of  data  restrictions. 

The  reason  for  this  is  rather  obvious.  The  purpose  of 
modelling,  in  the  first  place,  was  to  conserve  the  data 
base.  We  were  able  to  do  this  by  isolating  certain  re¬ 
lationships  which  permitted  the  consideration  of  a  multi¬ 
station  data  base.  Dew-point  spread  changes  do  not  lend 
themselves  to  such  considerations  since  the  local  in¬ 
fluences  which  determine  their  changes  cannot  be  separated 
out  from  the  rest  of  the  sample  by  simply  relying  on 
generalized  statistical  relationships.  Hence,  we  are 
caught  in  the  dilemma  of  being  able  to  employ  sophisticated 
methods  to  build  a  reliable  model  which  we  are  then  forced 
to  run  with  unreliable  data.  Another  possibility  is  to 
have  the  forecaster  provide  the  required  dew-point  in¬ 
formation  by  employing  other  techniques  and  materials 
available  to  him.  These  are  primarily  the  extrapolation 
of  immediately  past  trends  and  the  use  of  advective 
principles  and  deductive  informations  from  the  latest 
available  12-hour  prognostic  products.  Ideally,  this  pro¬ 
cedure  automatically  combines  the  best  of  climatology  with 
the  best  forecast  information  available  to  produce  a  happy 
marriage  of  the  two.  Hence,  we  set  forth  to  test  our 
procedures  in  the  operational  routines.  It  was  our  hope 
that  the  forecaster  could  weigh  the  scanty  climatic  evi¬ 
dence  provided  by  us  for  the  dew-point  spread  changes 
along  with  the  other  informations  available  to  him  to 
provide  the  required  dew-point  forecast  for  to  +  2.  Four 
Air  Weather  Service  stations  were  selected:  Norton  AFB, 
California,  Offutt  AFB,  Nebraska,  Andrews  AFB,  Maryland 
and  Warner  Robbins  AFB,  Georgia.  Each  station  was  provided 
statistics  such  as  those  shown  in  Figure  24a  and  b  to  test 
its  operational  routine.  A  digest  of  their  reaction  to 
the  technique  is:  A  knowledge  of  the  moisture  fields  and 
the  persistency  of  an  existing  ceiling  appears  to  be  use¬ 
ful  (in  post-mortem)  but,  unless  better  guidance  can  be 
provided  for  supplying  the  dew-point  spread  two  hours 
hence,  it  has  limited  operational  value. 

It  is  perhaps  also  fair  to  state  that  the  forecasters 
tended  to  view  the  materials  being  presented  more  as 
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Fig.  24a.  2-hour  conditional  probabilities 
for  initial  ceilings  less  than 
200  feet  (produced  as  in  figure  23) 
submitted  to  field  testing  at 
Andrews  AFB,  Maryland. 
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complete  solutions  than  as  climatic  aids.  Although  the 
field  evaluations  provided  less  causes  for  optimism 
than  we  had  hoped,  they  were  far  from  being  negative. 

In  fact,  we  were  much  encouraged  by  their  dedication, 
frankness  and  perception.  This  experience  merely  high¬ 
lighted  the  necessity  for  us  to  rework  our  technique 
to  circumvent  a  common  weakness  of  the  science-  in  the 
short-range  forecasting  problem,  little  information 
except  persistence  and  its  climatology  can  be  relied 
upon  in  making  a  forecast.  Therefore,  before  one  can 
incorporate  diagnostic  informations  at  tQ,  he  must  cir¬ 
cumvent  the  need  for  a  forecast  input  of  these  same 
criteria  at  t0  +  2.  This  we  have  proceeded  to  do.  Our 
means  of  doing  it  will  be  discussed  in  the  section  to 
follow . 

V.  A  MODEL  WHICH  INCORPORATES  DIAGNOSTIC  WIND  AND  MOISTURE 
INFORMATIONS  AT  to  WITHOUT  REQUIRING  THAT  THEY  BE 
PROVIDED  PROGNOSTICALLY  AT  tQ  +  2 

One  of  the  more  thought  provoking  statements  forth¬ 
coming  from  the  field  evaluation  was  supplied  by 
Captain  Strater  of  the  Andrew's  evaluation  team.  He 
questioned  at  length  the  wisdom  of  total  reliance  on 
dew-point  spread  changes  to  carry  the  meteorological 
history  of  ceiling  changes.  His  comment  was  much  deeper 
than  that  of  questioning  the  accuracy  of  the  temperature 
dew-point  values  as  they  are  presently  measured.  He 
particularly  emphasized  a  fact  that  we  had  overlooked. 

The  round-off  errors  introduced  in  reporting  temperatures 
and  dew  points  in  whole  degrees  automatically  provided 
"errors  "of  the  same  magnitude  as  those  found  in  the 
climatic  dew-point  statistics.  This  realization  has 
caused  us  to  rework  the  problem  by  deviating  from  the 
basic  concept  underlying  all  of  our  previous  models  as 
well  as  those  of  McCabe,  McAllister  and  Gringorten's  in 
a  very  fundamental  respect.  No  longer  can  we  compare 
the  uncond it ional  frequency  of  occurrence  value  at  one 
hour  to  that  of  a  subsequent  one  if  we  want  to  take  full 
advantage  of  the  observed  diagnostic  informations  avail¬ 
able  at  tQ.  Instead  we  must  demand  that  the  frequency  of 
occurrence  values  for  the  subsequent  hour  under  con¬ 
sideration  (in  this  case  tQ  +  2)  be  conditioned  on  the 
same  observed  moisture  and  wind  considerations  as  those 
which  were  used  for  the  initial  hour.  This  implies  a 
double  conditioning  process.  The  original  "unconditional" 
probabilities  at  tQ  +  2  now  become  conditioned  on  the 
existence  of  some  meteorological  parameter  at  tQ.  On 
those  occasions  when  ceilings  in  a  given  category  are  also 
being  observed,  these  occurrence  statistics  are  then  used 
to  determine  the  conditional  probability  that  the  event 
will  subsequently  recur  at  tQ  +  2. 
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The  procedure  is  as  follows:  The  unconditional 
frequency  of  ceiling  category  occurrences  for  each 
initial  hour,  to>  is  computed  exactly  as  before  when 
moisture  and  wind  stratification  procedures  were  in¬ 
voked.  However,  the  two  hour  later  frequency  of  oc¬ 
currence  statistics  are  based  only  on  that  fraction  of 
its  climatic  records  which  had  previously  met  the 
requirements  imposed  at  tp  except  for  the  existence  of 
a  particular  celling.  The  "unconditional”  probabilities 
as  they  are  now  being  computed  can  be  greater  than 
those  which  were  previously  considered  for  that  hour. 
This  fact  presumably  has  little  bearing  on  whether  or 
not  empirical  relationships,  with  the  same  universal 
tendencies  as  that  found  when  the  "purer"  statistical 
approach,  exist.  Since  this  supposition  demands  proof, 
we  are  proceeding  to  accumulate  it.  Preliminary  evi¬ 
dence  to  support  this  assumption  is  shown  in  Figure  25. 
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Fig.  25.  2-hour  statistically  recon¬ 
structed  conditional  probabilities 
based  on  a  "conditioned"  uncon¬ 
ditional  probability  at  tQ  +  2  and 
does  not  need  a  dew-point  spread 
forecast.  The  example  is  for 
Randolph  AFB.  Texas  under  calm  winds. 


-34- 


The  fundamental  advantage  to  the  technique  is,  of  course, 
the  elimination  of  a  requirement  to  forecast  the  dew¬ 
point  spread.  However,  the  model  previously  discussed 
in  Section  IV  can  still  offer  an  option  to  the  forecaster 
to  enter  such  information  should  he  feel  confident  in 
doing  so.  The  main  differences  between  the  models  of 
Section  IV  and  V  are  that  one  requires  an  explicit  dew¬ 
point  change  forecast  and  the  other  treats  this  parameter 
implicitly.  Since  the  observed  information  all  comes 
from  the  same  data  base,  Figure  23  and  Figure  25  should  be 
compatible  providing  the  "effective"  dew-point  spread 
changes  were,  in  fact,  known  for  insertion  into  Figure  23. 
Hence,  a  comparison  of  the  respective  modelling  results 
provides  a  means  of  inferring  the  effective  climatic 
dew-point  spread  changes. 

We  plan  to  proceed  in  this  direction  in  the  months  to 
come  and  prepare  climatic  presentations  for  reintroduction 
to  the  field  by  the  onset  of  the  1973-1974  winter  season. 

Our  plans  are  to  provide  the  climatic  forecast  as  a 
standard  and  let  the  field  forecaster  trend  his  other  rele¬ 
vant  materials  as  he  may  or  may  not  wish  to  do.  There  is 
little  doubt  in  our  minds  but  that  the  ultimate  procedure 
lies  in  some  such  man-climatology  mix.  The  primary  contri¬ 
bution  of  the  model  in  this  section  then  is  to  establish 
the  best  that  climatology  alone  has  to  offer  without  re¬ 
quiring  a  forecast  input  and  present  it  in  a  form  that  man 
can  provide  modifications  to  as  he  deems  fit  without  entirely 
abandoning  the  worth  of  climatology.  In  fact,  the  better 
the  man  can  provide  inputs  into  the  model  the  more  relevant 
the  climatic  information  becomes  to  the  specifics  of  the 
problem  at  hand.  The  need  for  this  mix  will  loom  in¬ 
creasingly  important  for  the  forecasting  of  rapidly  de¬ 
teriorating  ceiling  conditions  since  climatology  alone 
would  not  appear  to  suffice. 

We  have  not  abandoned  our  obligation  to  devise 
climatic  information  to  help  man  enter  these  forecast  in¬ 
puts.  Our  preliminary  investigations,  on  how  best  to  fore¬ 
cast  the  moisture  parameter,  show  that  the  problem  needs  to 
be  treated  in  two  parts.  The  temperature  component  of  the 
combination  conforms  quite  closely  to  the  types  of  climatic 
considerations  we  have  discussed  in  this  report  since  it  is 
strongly  affected  by  solar  insolations.  The  dew-point 
chtinge  component,  on  the  other  hand,  is  somewhat  less  de¬ 
pendent  upon  these  effects.  Here  past  weather,  its  trends, 
and  its  advecticn  from  local  source  regions  appear  to  play 
a  much  more  predominant  role  than  the  diurnal  climatolog¬ 
ical  amplitude  fluctuations  of  the  dew-point  curves. 

Attempts  at  forecasting  the  dew  point  independently  from 
the  temperature  should  automatically  lead  us  directly  into 
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one  of  the  main  goals  of  the  proposed  carry-on  re¬ 
search—  the  incorporation  of  synoptic  considerations 
common  to  given  regions  into  models  which  generate 
recurrence  probabilities  for  specific  point  locations 
within  that  complex. 

A  more  detailed  discussion  of  much  of  the  ma¬ 
terial  in  this  report  is  contained  in  an  article  pub¬ 
lished  in  the  Journal  of  Applied  Meteorology  (Martin, 
1972) . 

VI.  DESIGNING  CLIMATIC  INFORMATION  TO  SUPPORT  PARADROP 
OPERATIONS 

This  research  is  based  upon  the  following  assump¬ 
tions:  1)  Paradrops  are  released  at  heights  approxi¬ 

mating  1,000  feet  and  fall  freely  through  the  air  for 
some  150  feet  before  the  chute  opens,  2)  The  position 
of  the  aircraft  is  accurately  known,  3)  The  flight 
level  wind  direction  and  wind  speed  are  known  for  the 
release  point,  and  4)  The  surface  conditions  are  avail¬ 
able  either  through  remote  sensing  devices  or  from 
radio  contacts  with  a  ground  base  observer.  The  mete¬ 
orological  problem  is  to  determine  the  integrated  wind 
for  that  portion  of  the  trajectory  when  the  chute  is 
open.  A  guiding  criteria  is  for  the  wind  information 
to  be  sufficiently  accurate  that  it  leads  to  some  high 
percentage,  say  95%,  of  the  troops  landing  within  a 
prescribed  circular  drop  zone. 

The  vertical  structure  of  the  atmospheric  layer 
of  interest  covers  all  of  the  so-called  constant  flux 
layer  (some  50'  above  the  ground)  and  portions  of  the 
Ekman  layer  (a  theoretical  layer  which  extends  from 
the  top  of  the  constant  flux  layer  to  the  lower  extrem¬ 
ities  of  free  atmosphere). 

Theoretical  treatments  of  the  lower  layer  often 
carry  the  assumptions  that  the  earth's  topography  is 
relatively  level  and  quasi-homogeneous  and  that  the 
mean  pressure  gradient  is  invariant  with  height.  These 
conditions  signify  plane-parallel  motions  and  a 
constant-wind  direction  within  the  surface  layer.  The 
additional  assumption  of  quasi-steady  state  conditions 
restricts  the  surface  layer  to  one  of  constant  mean 
vertical  momentum  and  heat  fluxes.  According  to  the 
Monin-Obukhov  similarity  theory  (1954),  the  wind  speed, 
u,  and  temperature,  T,  profiles  in  this  layer  can  be 
expressed  by  the  dimensionless  universal  functions  f 
and  f T  :  U 
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u(z2>  -  =  3 r  [fu(  L1) '  Mr1)] 

T(z2>  -  T(zl>  =  T*  [fT  ( )-  fT(  IP ;'] 


where  z^ ,  z 2  are  arbitrary  levels,  k  is  the  von  Karman's 
constant,  u*  is  the  frictional  velocity,  T*  and  L  are 
the  characteristic  temperature  and  length  scales  defined 
as 


T* 


Here  q  is  the  vertical  turbulent  heat  flux.  fD0  and  T0 
are  "standard  values"  of  air  density  and  temperature  in 
the  surface  layer.  Cp  is  the  specific  heat  for  dry  air 
at  constant  pressure  and  g  is  the  gravitational  acceler¬ 
ation. 

Relatively  few  investigations  have  been  made  on 
the  wind  characteristics  for  the  Ekman  layer.  The  best 
known  model  is  due  to  Taylor  (1916)  which  is  based  on 
the  assumption  of  steady  motion,  no  advection,  neutral 
stability,  no  thermal  wind,  and  constant  eddy-viscosity 
coeff icient s . 

A  widely  accepted  empirical  relationship  which  re¬ 
lates  the  wind  speeds  at  different  heights  is  the  so- 
called  power  law, 


This  law  states  that  the  wind  speeds  at  different 
levels  can  be  fitted  by  a  straight  line  with  a  slope  of 
1/p  on  a  log-log  diagram.  Unfortunately  the  exponent, 
p,  varies  with  stability  conditions  and  geography 
(Davenport,  1960),  and  the  theoretical  wind  spirals  be¬ 
low  3,000  feet  often  display  little  or  no  resemblance 
to  the  observed  evidences  (Blackadar,  1960). 

Our  approach  to  the  problem  of  estimating  the  inte¬ 
grated  wind  for  the  paradrop  problem  at  hand  was  to 
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determine  statistical  relationships  between  the  inte¬ 
grated  wind  and  various  combinations  of  flight  level 
and  surface  wind  informations.  These  statistics  were 
then  used  to  provide  estimates  of  the  integrated  wind 
through  the  layer  intervening  between  the  two  levels. 

Two  data  sets  were  obtained  from  instrumented 
towers.  One  of  these  was  Cedar  Hill,  Texas,  and  the 
other  was  in  Oklahoma  City,  Oklahoma.  The  Cedar  Hill 
Tower  is  located  approximately  20  miles  SSW  of  Dallas, 
Texas,  at  an  elevation  of  850  feet  above  sea  level. 

The  terrain  elevation  decreases  by  about  300  feet  over 
a  distance  of  3  to  5  miles  to  the  west  and  north  of 
the  site.  The  topography  is  essentially  flat  to  the 
south  and  east  (see  Gerhardt  et  al,  1962). 

The  instrumentation  consists  of  Bendix-Friez 
Aerovanes  (t  3°  in  direction)  and  copper-constantan 
thermocouples  (  \ error!  <  1°F)  located  at  12  levels: 

30,  70,  150,  300,  450,  600,  750,  900,  1,050,  1,200, 

1,300  and  1,420  feet.  These  sensors  are  mounted  on 
booms  extending  12  feet  from  the  tower  structure.  A 
continuous  record  of  digital  recorded  10  minute-averaged 
wind  and  temperature  observations  on  magnetic  tape  for 
December,  1960  to  December,  1962  was  used  in  this  study. 

The  Oklahoma  City  data  were  taken  on  the  instru¬ 
mented  WKY-TV  tower.  The  terrain  about  this  site  is 
gently  rolling  with  an  average  height  of  1,148  feet 
above  sea  level.  A  Bendix  Aerovan  wind  sensor  and 
thermistor  (±  1°C)  were  used  to  take  the  observations. 

The  instruments  are  mounted  at  six  levels:  146,  296,  581. 
874,  1,166  and  1,458  feet.  A  surface  station  250  feet 
away  provided  concomitant  wind  data  averaged  over  a 
period  of  five  minutes  every  hour  on  the  hour  at  the  23- 
foot  level.  The  wind  data  from  this  tower  which  was 
used  in  our  study  extends  from  May,  1966  to  May,  1967. 

The  temperature  registrations  were  somewhat  of  lesser 
duration  extending  from  December,  1966  to  May,  1967. 

For  further  details  with  respect  to  the  data  and  its 
accuracy  see  Carter,  1970. 

The  procedure  which  we  used  for  calculating  the  in¬ 
tegrated  winds  from  observed  data  is  as  follows.  The 
u  and  v  components  of  the  integrated  wind  were  obtained 
as  a  weighted  average  of  the  u  and  v  components  of  the 
observed  winds  at  each  of  the  reporting  levels.  For 
example,  the  observed  wind  at  600  feet  was  assumed  to 
be  representative  of  the  525  to  675  foot  layer.  The 
weighting  factor  for  this  wind  is  150  feet.  Since  the 
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thickness  of  the  given  layer  is  a  function  of  the 
spacing  between  the  instrumented  levels,  the  weighting 
factors  vary  throughout  the  vertical. 

The  winds  at  the  6-foot  level  were  found  by  ex¬ 
trapolating  downward  from  the  adjacent  highest  level 
using  the  power  law 


Singer  et  al  (1970)  estimated  the  exponent  p  to 
have  a  value  of  0.2.  We  used  his  value  since  the  rela¬ 
tively  small  height  increments  through  which  these  ex¬ 
trapolations  were  made  did  not  appear  to  warrant  the 
introduction  of  diurnally  or  seasonally  varying  ex¬ 
ponents  . 

Some  diurnal  and  seasonal  characteristics  of  the 
integrated  wind  speeds  at  each  of  the  contributing  levels 
in  the  Cedar  Hill  data  are  found  in  figures  26  through  29. 


L?T  22  0  *  *  i  •  10  12  in  16  II  20  22 

Tig,  26.  Mean  wind  speed  distributions  at  various 
levels  for  Cedar  Hill  during  January. 
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Fig,  37.  Mean  wind  speed  distributions  at  various 
level*  for  Cedar  Hill  during  April. 


Fig.  31.  Mesn  wind  speed  distributions  at  various 
level*  for  Cedar  Hill  during  July. 
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Tig.  21 .  Mean  wind  tp««d  distributions  at  various 
lavals  for  Cedar  Hill  during  October. 


The  averaged  winds  for  each  level  were  computed  from 
the  data  of  that  10-minute  time  interval  which  ended 
1  minute  before  the  hour.  The  data  in  these  figures 
show  the  integrated  wind  to  correlate  closely  to  a 
mid-layer  wind  and  that  this  wind  is  not  an  algebraic 
average  between  the  flight  level  and  surface  winds. 

Note  that  the  respective  wind  profiles  at  levels  above 
450  feet  show  more  homogeneity  among  each  other  than 
do  the  winds  at  levels  below  these  heights  and  that 
the  integrated  wind  profile  exhibits  its  best  phase  and 
amplitude  relationship  with  the  higher  grouping.  The 
inference  is  rather  apparent .  The  most  single  per¬ 
tinent  piece  of  meteorological  information  which  could 
be  entered  into  the  parachute-trajectory  equations 
would  be  a  wind  observation  midway  between  the  paradrop 
release  point  and  the  ground.  Thus,  it  is  recommended 
that  such  a  pass  be  taken  whenever  it  is  feasible  to 
do  so.  When  and  if  such  procedures  are  prohibitive, 
the  most  valuable  information  is  the  flight-level  wind 
at  the  release  point. 
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The  techniques  that  we  shall  discuss  apply  only  when 
the  first  alternative  (mid-level  observation)  is  not  feas¬ 
ible.  Should  such  mid-level  observations  be  available,  we 
recommend  that  the  empirical  methods  which  we  have  devised 
for  estimating  the  integrated  wind  be  supplanted  by  them. 

Many  of  our  attempts  either  proved  fruitless  or  became 
too  involved  for  eventual  operational  application.  They 
will  not  be  included  since  our  objective  is  to  seek  posi¬ 
tive  results  without  becoming  overly  involved  in  tidying  up 
our  failures  for  historical  purposes.  Hence,  when  any  par¬ 
ticular  line  of  attack  showed  overwhelming  evidence  of  not 
being  able  to  meet  certain  criteria  it  was  immediately  dis¬ 
continued  irrespective  of  the  efforts  we  had  spent  upon  it. 
These  criteria  are:  Will  a  continuance  of  the  approach  pro¬ 
duce  a  better  climatic  product  than  those  already  at  hand, 
and  will  the  end  results  be  sufficiently  straightforward 
and  readily  attainable  that  they  can  be  applied  within  the 
constraints  imposed  by  real-time  operations? 

The  simplest  approach  we  have  found  which  meets  these 
aforementioned  criteria  relies  solely  on  the  latest  ob¬ 
served  wind  at  flight  level.  It  is  presented  in  Figures  30 
through  37  for  the  months  of  January  and  July.  These  were 
obtained  by  partitioning  the  Cedar  Hill  winds  for  each  month 
at  the  assumed  flight  level  of  1,050  feet  into  increments  of 
5  mph  speed  and  450  direction. 
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Fij.  JO.  Deviation  vectors  for  calculating  the  inte¬ 
grated  wind  at  Cedar  Hill,  Texas  for 
2100-02S9  during  January. 
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Flj.  31.  Deviation  vectors  for  calculating  the 
integrated  wind  at  Cedar  Hill,  Texas 
for  0300-08 S9  during  January. 


rig.  37.  Deviation  vector*  for  calculating  the 
integrated  wind  at  Cedar  Hill,  Texas 
for  0900-1459  during  January. 
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770° 


rig.  i).  Deviation  vector#  for  calculating  the 
integrated  wind  at  Cedar  Mill,  Texa# 
for  1500-2059  during  January. 


270° 


rig.  3% 


Deviation  vectors  for  calculating  the 
integrated  wind  at  Cedar  Hill,  Texa# 
for  2100-0253  during  July. 
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Tig.  3$.  Deviation  vectors  for  calculating  the 
Integrated  wind  at  Cedar  Hill,  Texas 
for  0300-08S9  during  July. 


110° 


rig.  36.  Deviation  vectors  for  calculating  the 
integrated  wind  at  Cedar  Mill,  Texas 
for  0900-X1S9  during  July. 
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Fig.  37.  Deviation  vector*  for  calculating  tha 
integrated  wind  at  Cedar  Hill,  Texa* 
for  1500-2059  during  July. 

The  observations  at  each  respective  lower  level 
were  individually  processed  to  obtain  the  integrated 
wind  which  pertained  to  each  respective  flight  level 
wind  observation.  These  integrated  winds  were  then 
vector-averaged  to  obtain  a  mean  for  the  particular 
subset  that  it  represents.  These  presentations  permit 
a  graphical  estimation  of  the  integrated  wind  from  a 
knowledge  of  the  flight  level  wind.  For  example, 
assume  that  the  observed  flight  level  wind  at  0630  local 
solar  time  in  January  was  225°  at  22.5  mph.  The  inte¬ 
grated  wind  is  found  by  drawing  a  vector  from  the  center 
of  the  display  to  the  head  of  the  arrow  eminating  from 
the  2250  and  22.5  mph  location.  A  value  of  some  220° 
and  18  mph  is  readily  attained.  Identical  procedures 
applied  to  other  wind  direction  for  the  22.5  mph  wind 
speed  value  produces  integrated  winds  of  relatively 
the  same  magnitude  (18  mph).  However,  the  directional 
relationships  between  the  upper  level  and  integrated 
winds  are  found  to  be  somewhat  more  dependent  on  the 
initial  wind  direction.  Referring  to  Figure  26,  a 
value  slightly  greater  than  17  mph  is  given.  Thus, 
the  displays  in  Figs.  30  through  37  add  another  im¬ 
portant  dimension  (initial  wind  direction)  to  that  given 
by  Figures  26  through  29. 

The  reliability  of  the  "deviation  vectors"  on  these 
graphs  is  obviously  a  function  of  the  sample  size. 
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In  our  research  we  discarded  all  data  points  having 
less  than  5  cases  for  constituting  the  average.  For 
example,  the  wind  at  Cedar  Hill  has  a  predominantly 
southerly  component  during  July.  The  rarely  occurring 
northerly  components  account  for  the  general  blank  area 
about  the  360°  sector. 

We  recommend  this  simple  climatic  presentation  for 
operational  implementation  whenever  the  winds  at  flight 
level  constitute  the  only  meteorological  input  into  the 
parachute-traj ectory  computation . 

The  best  approaches  which  we  found  to  apply  when 
the  6-foot  wind  direction  and  speed  were  also  known  of 
the  integrated  wind  will  be  presented  at  this  time. 

The  first  one  to  be  discussed  involves  regression  equa¬ 
tions  of  the  form 

VI  5  ao  t  al  v6  *  a2  v1050 

where  vg  and  v^Qg0  are  the  wind  speeds  at  6  feet  and 
flight  level  respectively,  the  vj  is  the  integrated 
wind  for  that  portion  of  the  intervening  layer  where  the 
chute  is  assumed  to  be  open  (825-15*).  The  respective 
equations  for  Cedar  Hill  are  shown  in  Figure  38.  Here 
for  convenience,  the  standard  errors  of  estimate  are 
labeled  as  root  mean  square  errors  (RMSE).  They  are 
found  to  be  the  smallest  during  the  time  period  09-15 
for  all  seasons  since  this  interval  is  conducive  to  a 
maximum  of  vertical  coupling  and  momentum  exchange  be¬ 
tween  the  two  levels. 

An  alternate  means  for  estimating  the  integrated 
wind  from  a  knowledge  of  those  at  flight  level  and  the 
ground  is  shown  in  Figures  39  through  42.  The  main  ad¬ 
vantages  of  this  technique  over  the  regression  equation 
approach  are  the  relaxation  of  some  of  the  linearity 
constraints  and  a  presentation  of  the  relationships  as 
a  field  variable  more  conducive  to  pictorial  assessments. 

To  illustrate  this  point  select  isolines  from  each  of 
these  four  figures  were  transposed  onto  one  graph 
(see  Figure  43).  Here  it  is  seen  that  the  specific  com¬ 
binations  of  flight  level  and  surface  winds  which  produce 
estimated  integrated  winds  on  the  order  of  5  mph  are 
relatively  invariant  throughout  the  day  for  this  location 
and  month.  A  similar  statement  would  pertain  to  higher 
values  equal  to  or  greater  than  25  mph.  The  fact  that 
the  intermediate  values  (as  represented  by  the  15  mph 
isolines)  show  the  most  variance  immediately  serves  to 
alert  the  user  to  those  cases  when  a  climatic  approach  to 
the  problem  is  the  least  likely  to  give  acceptable  results. 
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ft*.  1*.  lnt«|nt«4  Wind  S»««d  (1S-»?I  Tt.)  in  H*H 
Coder  fill  {October  J100-02S»> 


rii.  *0 ,  Inteftreted  Wind  5p««d  (U-I2S  T t.)  in  HPK 
Ce<J*r  Hill  (October  3130  -  0*S3) 
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Fig.  41.  Interpreted  Wind  Speed  (15-1??  rt.)  in  KFK 
Coder  HiU  (  October  OIOO-IOI) 


rif.  4J.  Int tf.MtH  Wind  Speed  ( rt.>  in  wr 
Ceder  Hill  (  October  1700  -  70S!) 


S|4e«,d  (In  Hill!  *t  *.  r-nt 
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rir.  «1,  Int**r*t*tf  Wind  Pp**4  Ft.)  In  HPH 

C«d«r  Hill  <  October)  for  Various 
•-Hourly  Tl»*  Periods 


In  the  final  analysis,  the  value  of  each  of  the 
aforementioned  climatic  aids  is  dependent  upon  the 
"universality"  of  the  relationships  which  they  contain 
since  it  is  highly  unlikely  that  separate  dependent  data 
displays  will  ever  become  available  for  every  conceivable 
release  point  throughout  the  world.  Thus,  it  is  neces¬ 
sary  that  these  aids  be  tested  on  independent  data  to 
determine  the  nature  and  extent  of  the  given  relation¬ 
ships.  We  applied  the  regression  coefficients  which  were 
derived  from  the  Cedar  Hill  Tower  data  to  the  flight 
level  and  surface  wind  data  of  Oklahoma  City  to  produce  an 
estimated  integrated  wind,  Vj.  The  results  are  shown  in 
Figure  44.  RMSE,  the  standard  error  of  estimate,  pertains 
when  tfte  data  for  Oklahoma  was  independent  to  the  re¬ 
gression  coefficients.  The  term  RMSE*  pertains  when  the 
data  from  Oklahoma  City  were  used  to  derive  another  set 
of  coefficients  (not  shown)  to  obtain  a  dependent  esti¬ 
mate  of  the  integrated  wind  v^  The  term  RMSE/vj 
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3.367 

V10S0 

11.271 

1.195 

11.051 

11.311 

*1 

10.5SS 

1.111 

1.125 

10.115 

*1 

13.059 

6.155 

10,012 

1}.706 

RMSE 

2.510 

1.306 

1,266 

2.506 

RMSE/Vj 

0.231 

0.156 

0.1)6 

6.226 

RMSE* 

1.111 

0.712 

0.651 

1.156 

Fift.  11.  Independent  and  deomdent  t*«ta 
of  the  rerremion  equation 
VI  '  *o  *  *ly10$« 
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tends  to  normalize  the  independent  data  error  according 
to  the  magnitude  of  the  integrated  wind.  It  is  logical 
that  the  RMSE  for  independent  data  would  exceed  the 
RMSE*  of  the  dependent  set.  Nevertheless,  in  many  cases, 
particularly  for  the  09-15  and  15-21  time  periods  for 
each  season,  the  comparative  magnitudes  of  the  two 
terms  are  similar. 

Presently  the  integrated  winds  are  calculated 
operationally  by  averaging  the  surface  and  flight  level 
wind  (see  AWS  special  study,  1971).  Figure  45  compares 
the  standard  errors  of  estimates  for  the  aforementioned 
independent  data  verification  for  the  Cedar  Hill  re¬ 
gression  coefficients  on  the  Oklahoma  City  data  versus 
that  which  would  have  been  attained  had  the  simple 
averaging  scheme  been  employed  in  its  stead.  These  sta¬ 
tistics  indicate  that  the  regression  technique  is  rela¬ 
tively  more  reliable  during  the  hours  of  09-2100  for 
all  seasons.  Hence,  if  paradrop  operations  are  con¬ 
fined  to  these  hours,  the  technique  we  have  devised 
offers  possible  improvements  over  the  operational  pro¬ 
cedures  currently  in  existence  during  those  times  when 
active  mixing  processes  are  in  being.  The  fact  that 


Local  Standard  Tine 


January 

03-09 

09-15 

1S-21 

21-03 

a. 

k.lkS 

2.386 

3.62k 

3.783 

b. 

3.068 

2 . 095  • 

2.122 

2.2k3 

April 

a. 

2.987 

3.586 

3.k09 

3.189 

b. 

k.105 

1.935 

1.6S0 

3.2k8 

July 

a. 

1.133 

1.08  9 

1.380 

1.339 

b. 

October 

1.291 

0.638 

0.961 

2.362 

a« 

1.933 

1.S96 

1.900 

2.076 

b. 

2  •  SkO 

1.306 

1.266 

2.S09 

Note: "a. 

V 

"represents 
I  *  1/2  <V« 

RMSE  due  to  VT 
+  V10S0> 

estimation 

using 

"b. "represents  RMSE  due  to  V_  estimation  when  th- 
regmesion  coefficients  developed  for  the 
Cedar  Hill  Tower  data  wnre  used  to  make  ir.tee.r^ted 
wind  estimates  on  the  data  base  for  Oklahoma  City. 

rig.  kl.  Comparison  of  standard  errors  of  estimates  in 
obtaining  Vj,  Oklahoma  City  Tower. 
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the  "nighttime"  hours  verify  somewhat  poorer  on  the  in¬ 
dependent  data  set  than  the  daytime  hours  suggests 
that  stability  considerations  be  introduced  more  directly 
into  the  method.  This  parameter  was,  indeed,  found  to 
provide  a  decided  improvement  in  the  accuracy  of  esti¬ 
mation  when  the  observed  temperature  profile  throughout 
the  layer  was  known.  However,  a  linear  representation 
of  the  profile,  found  by  simply  considering  the  tempera¬ 
tures  at  flight  level  and  the  surface,  did  not  suffice. 
Hence,  the  presentations  were  altered  to  better  in¬ 
corporate  the  diurnal  stability  trends  implicitly  by 
increasing  the  diurnal  resolution.  The  need  for  this 
is  clearly  seen  from  an  examination  of  the  statistics 
presented  in  Figure  46. 


Hour 

Number  of  Data 

0-2 

257 

2-4 

273 

4-6 

286 

6-6 

278 

8-10 

235 

10-1? 

260 
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14-16 
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16-10 

233 

18-20 

297 

20-22 

30? 

22-24 

28S 
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3.0 

3.4 

3.3 

2.0 

1.4 

-0.3 

-1.0 

-0.5 

0.2 

1.4 

3.0 

3.1 

al 

0.54 

0.53 

0.53 

0.50 

0.41 

0.27 

0.31 

0.46 

0.46 

0.40 

0.45 

0.47 

*2 

0.48 

0.41 

0.47 

0.74 

0.98 

1.30 

1.30 

0.98 

0.98 

1.22 

0.83 

0.66 

RMS 

2.21 

2.26 
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2.47 

2.66 

1.13 

0.99 

0.89 

1.75 

2.64 

2.17 

2.26 

Mean 

of  V10S0 

20.1 

20.2 

19.2 

17.2 

16.3 

13.3 

11.6 

11.2 

13.9 

17.0 

16.3 

19.4 

Std .  dev. 

of  v1050 
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7.4 
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8.9 

8.6 

8.5 

7.9 

0.1 

8.0 

7.3 

Mean 

of  v6 

6.3 

5.0 

5.8 

5.3 

6.2 

6.0 

5.4 

5.2 

5.2 

4.8 

4.9 

5.2 

Std.  dev. 

of  v6 

2.5 

2.6 

2.5 

2.  6 

3.4 

3.7 

3.9 

3.7 

3.0 

2.1 

2.1 

2.3 

Mean 

of  Vj 

16.9 

16.5 

16.3 

14.7 

14.1 

11.1 

9.7 

9.7 

11.8 

14.1 

15.3 

15.7 

Std.  dev. 

of  V 

5.4 

5.4 

5.5 

5.6 

6.7 

6  .9 

7.6 

7.5 

6.5 

5.9 

5.5 

5.0 

Fig.  46.  Regression  coefficients  for  Cedar  Hill 
during  January  using  th«  equation 

vi  *  ao  ♦  ai  vioso  *  *2  ve 


The  diurnal  variation  of  a2  is  found  to  be  more  than 
300%  and  that  of  a^  some  200%.  The  a]_  coefficient  shows 
a  uniform  diurnal  bias  with  a  minimum  near  mid-day 
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similar  to  that  evidenced  by  the  integrated  wind  profile 
(see  Figures  26  through  29).  The  a-2  coefficient  on  the 
other  hand  shows  a  lesser  diurnal  regularity.  Thus, 
a  six-hourly  stratification  of  the  data  might  be  warranted 
in  the  derivation  of  a  regression  equation  if  upper  level 
winds  alone  were  to  be  considered.  However,  this  time 
interval  is  much  too  long  for  partitioning  the  data  when 
surface  wind  considerations  are  also  introduced.  Since 
hourly  or  bi-hourly  stratifications  would  be  rather  un¬ 
wieldy  and  unreliable,  due  to  a  shrinking  of  the  data 
base,  we  propose  that  other  means  than  the  regression 
equation  be  sought  for  introducing  the  surface  wind  in¬ 
formation  into  the  integrated  wind  estimates.  Some  of 
these  possibilities  are  in  the  checking-out  stage. 

Attempts  were  also  made  to  devise  schemes  for  es¬ 
timating  the  flight  wind  at  a  given  point  location 
from  a  knowledge  of  the  observed  winds  at  nearby  lo¬ 
cations.  Rawin  data  from  Fort  Worth,  Oklahoma  City, 

Little  Rock,  Shreveport,  Lake  Charles,  San  Antonio, 

Midland  and  Amarillo  were  used  in  various  combinations 
to  estimate  the  1,300  foot  winds  for  Cedar  Hill.  For 
example,  one  experiment  considered  the  data  of  all  eight 
stations,  another  only  that  of  a  few,  etc.  The  best 
set  of  correlation  statistics  were  found  when  only  wind 
information  for  the  geographically  most  proximate  station, 
Fort  Worth,  was  considered.  A  21-dav  test  during  January 
showed  this  regression  equation  to  be  y  =  x  -  2  where 
y  is  the  1,100  foot  wind  at  Forth  Worth  and  x  is  the 
1,300  foot  wind  at  Cedar  Hill.  The  RMSE  of  this  estima¬ 
tion  was  less  than  4  mph. 

If  the  wind  distribution  is  assumed  to  be  circular 
normal,  then  the  u  and  v  components  are  independently 
distributed  and  the  standard  deviation  of  wind  speed  is 
related  to  the  variance  of  the  components  by 


The  formula 


describes  the  circle  of  radius  or  that  contains  the  tip 
of  the  wind  vector  with  probability,  P. 

Figure  47  represents  the  climatic  mean  wind  vector 
for  v^Q5o  (at  Cedar  Hill  for  09-15  October)  that  has  a 
speed  range  10-15  mph  and  a  direction  range  157 . 5°-202 . 5°. 


ITO*“ 


fi*.  *7. 


Probability  distribution  of  th«  10^0  foot  wind 
(b*«*d  upon  JJJ  et«««)  for  C*d*r  Rill  durinr  th* 
hour*  OSOO-liOO  in  October  wh*n  thwt  wind  was 
•b**rw*d  b»tw*«n  11?. S  nnd  707, $  dc#r««>*  with 
•»**d*  b«tw*nn  13-15  oph,  Th*  elto«ti«  •< *n 
(find  v*«tor  *nd  probability  eirxi**  *r«  shown. 


A  mean  vector  wind  of  184°,  11.3  mph  was  found  from  a 
sample  of  222  cases.  The  number  on  the  circles  denote 
the  percentage  probability  that  the  tip  of  the  wind 
vector  would  fall  on  or  within  that  circle. 

Another  important  consideration  in  this  work  in¬ 
volves  the  time  variability  of  the  wind  information. 

For  example,  if  and  V2  denote  two  random  variables 
for  an  initial  time  and  a  subsequent  time  respectively, 
their  interdependence  may  be  determined  from  their  co- 
variance  tensor 


V 


I 

li 


Here  i  and  j  represent  the  components  of  the  vectors 
and  the  prime  denotes  the  deviation  from  the  sample  mean. 
The  problem  is  to  best  estimate  V2  from  a  given  By 

assuming  the  regression  is  linear  and  applying  the  least 
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square  criterion,  this  best  estimate  is  expressed  as 

the  scalar  V21  =  B^j  V^j  providing  the  autogression _ 2 

tensor  B^j  satisfies  the  requirement  that(V2i  -  Bij Vij ) 
is  a  minimum. 

The  elements  of  B^  are  obtained  by  setting  the 
derivative  of  this  expression  with  respect  to  B^j  equal 
to  zero.  The  first  element  of  the  autogression  tensor  is 
found  to  be 


Other  elements  are  similarly  derived. 

The  curves  in  Figures  48  and  49  were  derived  from 
two  years  of  wind  data  at  Cedar  Hill.  Ten  minute  averages 
of  u’s  and  v's  were  used.  The  standard  error  of  estimate 
(for  convenience  labeled  as  RMSE),  the  mean  phase  error 
and  the  mean  relative  errors  in  establishing  subsequent 
wind  speeds  at  1050  km  from  a  knowledge  of  the  initial 
ones,  at  0100  and  1300  in  January  are  presented  as  a 


fig.  M.  Autoragr^aalon  for  10$0  foot  winds  at 
Cadar  Hill,  Taxaa  with  Initial  tUa 
0100  LST  during  January. 
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T i«.  M.  Autor**r««* ion  of  10$Q  foot  wind* 
«t  Ce<*«r  HU1,  T***f  with  initial 
tU«  1100  LST  during  J«nu«ry. 


function  of  lag  time.  Here  the  mean  phase  error  (MPE) 
is  defined  as  the  average  of  the  absolute  directional 
differences  between  the  observed  wind  and  the  estimated 
winds  for  a  particular  lag  time.  The  mean  relative 
error  is  the  average  ratio  of  the  errors  in  speed  esti¬ 
mations  to  the  observed  wind  speed.  Each  set  of  errors 
is  seen  to  grow  most  rapidly  during  the  first  10  to 
20  minutes  after  observation. 

Figures  50  and  51  show  a  similar  set  of  statistics 
for  the  integrated  wind.  This  wind  is  seen  to  exhibit 
a  somewhat  better  set  of  statistics  (errorwise)  than  the 
flight  level  wind.  These  illustrations  simply  confirm 
the  fact  that  the  accuracy  of  wind  information  over  the 
drop  zone  perishes  rapidly  with  time.  Every  attempt 
should  be  made,  therefore,  to  confine  the  time  difference 
between  observation-  and  release-time  to  less  than  30 
minutes . 
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fig.  M  Autonr.niilan  tor  irt*rr«t»o  wind* 

€a4ar  Mill,  T*mi  with  initial  tlM  0100 

L.*T  4wrih«  Jtiwary. 


rig.  II.  Autorarrraaion  for  lntarratad  wm« 

at  Cadar  NUl,  ?a*aa  with  *nitl«l  ri»* 
1100  LST  Aurliui  January. 
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The  directional  deviation,  d,  is  defined  as  the 
angular  difference  between  the  flight  level  wind  and 
the  integrated  wind.  A  positive  (negative)  deviation 
indicates  that  the  flight  level  wind  veers  (backs) 
from  the  integrated.  If  the  idealized  Ekman-Taylor 
spiral  representation  were  to  pertain  for  the  wind  dis¬ 
tribution  between  the  flight  level  and  the  ground,  then 
"d"  would  be  a  positive  value.  Our  empirical  results 
indicate  that  a  first  approximation  to  the  integrated 
wind  direction  would  be  a  direction  of  10°  less  than 
the  observed  flight  level  wind.  A  pictoral  assessment 
of  those  initial  wind  directions  which  give  integrated 
winds  that  deviate  most  from  this  average  is  provided 
by  the  arrows  in  Figures  30  through  37.  A  statistical 
breakdown  of  the  directional  deviations  which  lie 
between  -5  and  +25°  for  different  hours  and  seasons 
for  the  Oklahoma  and  Cedar  Hill  data  is  given  in 
Figure  52. 

VII.  SUMMARY  AND  CONCLUSIONS 

This  research  has  demonstrated  that  improved  cli¬ 
matic  aids  can  be  provided  for  the  forecasting  of  low 
ceilings  and  restricted  visibilities.  It  was  found 
that  the  conditional  probabilities,  which  are  currently 
being  provided  to  the  field  forecaster,  are  misleading 
for  the  less  frequently  occurring  events.  This  diffi¬ 
culty  can  be  partially  overcome  by  subjecting  these 
probabilities  to  the  hour-and-month  continuity  con¬ 
siderations  discussed  in  this  report.  This  study  also 
shows  that  conditional  probabilities,  of  comparable 
quality  to  those  currently  being  produced  by  computer¬ 
tracing  methods,  can  be  provided  by  empirical  formula¬ 
tions  on  a  much  smaller  data  base  than  presently  needed. 

The  persistency  of  a  cloud  ceiling  is  a  function  of 
the  wind  and  moisture  conditions  occurring  in  associ¬ 
ation  with  it.  Methods  for  introducing  these  considera¬ 
tions  into  the  recurrence  statistics  were  demonstrated 
for  various  degrees  of  modelling  sophistications. 

The  integrated  winds,  characteristic  of  a  drop- 
zone  layer,  were  studied  as  functions  of  the  flight 
level,  the  surface  wind,  and  various  diurnal  stability 
processes.  Several  different  empirical  methods  were 
presented  for  estimating  the  integrated  winds. 

Our  future  efforts  will  adapt  these  procedures  to 
the  4-hour  time  period  and  apply  them  to  the  prediction 
of  visibility  as  well. 
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Each  of  the  2-  and  4-hour  methodologies  will  be 
tested  in  a  two  fold  manner.  One  will  use  the  entire 
data  records  of  several  independent  stations  as  a 
verfication.  The  other  will  introduce  the  new  sets  of 
universal  graphs  into  the  field  to  determine  their 
operational  acceptability. 
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Note:  (a)  Range  of  integration  for  V^:  CHT  (1)  15-825  ft., 

CHT  (2)  15-1125  ft.,  OKC  (1)  11.5-874  ft., 

OKC  (2)  11.5-  11GB  ft. 


Fig.  52 


The  percentage  of  directional  deviations 
(d)  which  lie  between  -5°  and  +25°  for 
different  hours  and  seasons  for  the 
Oklahoma  (OKC)  and  the  Cedar  Hill  (CHT)  data. 
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13.  ABSTRACT 

This  research  considers  the  devising  of  climatic  presentations 
for  usage  in  preparing  2-  and  4-hour  forecasts  of  low  ceilings  and 
restricted  visibilities  and  the  providing  of  integrated -wind  infor¬ 
mation  within  the  drop-zone  layer  in  support  of  paradrop  operations. 

A  hierarchy  of  empirical  relationships  for  estimating  the  re¬ 
currence  probabilities  under  a  number  of  contingencies  are  made 
possible  through  the  use  of  statistical  modelling  techniques  which 
act  to  conserve  the  data  base.  These  are:  1)  a  model  for  genera¬ 
ting  recurrence  probabilities  which  considers  only  the  latest  ob¬ 
served  ceiling  category,  2)  one  which  incorporates  surface  wind  in¬ 
formation,  3)  another  which  includes  moisture  considerations  and  4) 
one  which  introduces  observed  moisture  and  wind  information  diag¬ 
nostically  without  having  to  assume,  the  persistency  of  these  vari¬ 
ables,  or  otherwise  forecast  their  behavior,  throughout  the  fore¬ 
cast  interval. 

Climatic  presentations  for  estimating  integrated  winds  within 
the  drop-zone  layer,  from  a  knowledge  of  the  winds  at  flight  level, 
are  discussed  apart  from  those  which  pertain  when  surface  winds  are 
similarly  available. 
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